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1 Introduction 



In spacetime geometry scenarios of general relativistic astrophysics and cosmology it often happens that 
there exists a preferred timelike vector field. For example, one might have a source corresponding to a stress- 
energy-momentum tensor with a timelike eigendirection, perfect fluids being perhaps the most familiar 
example. Other cases arise when symmetries give rise to a preferred timelike direction. Spatially self-similar 
and spatially homogeneous models are examples of this situation. If there exists a preferred timelike vector 
field on the spacetime manifold then the associated unit timclikc vector field u (such that u'' = — 1) 
determines tensors U^i, and h'^i,, respectively projecting parallel and orthogonal to u, according to 



It follows that 



(1.1) 

U%UP, = , U'^.u'^ = , [/"^ = 1, , . 

hi'p hP, = hi",, hi^^u" = 0, hi"^ = 3. ^ ■ ' 

Due to the structure imposed on the spacetime manifold by the existence of u, a 1+ 3 tensor decomposition of 
all geometrical objects of physical interest can be made with the help of the projection tensors If^^, and h^^,. 
For example, it is standard to decompose the covariant derivative Vj^m^ into its irreducible parts according 
to [1] 

V^Mi. = -UfiU„ + (Tfiu + 5 hfiu - OJflu , (1-3) 

where u is symmetric and tracefree, u!^,ii is antisymmetric, and — iifj^u^ — Ofinvy — LOfmvy. The 
kinematical fields associated with the timelike congruence u are defined by 



(1.4) 



where ii'^ denotes the acceleration vector; the (volume) rate of expansion scalar; (t^,, is the rate of shear 
tensor, with magnitude 

(7^ := i (7^. a'"' , (1.5) 
and uj^i, is the vorticity tensor. It is convenient to define a vorticity vector 

:= i r?'"^"'^ ^ LOf,, = ri,,p, lj" , (1.6) 

where rj^-"'"' is the totally antisymmetric permutation tensor. The sign convention jf^'^^ = 1/^/3^^ ??oi23 = 
— (where g denotes the determinant of the spacetime metric tensor, g^^y, with Lorentzian signature), 

will be employed throughout, following the conventions established in Ref. [2]^. The magnitude u) of the 
vorticity is defined by 

LO^ ■.= UJ^UJ^ = ^10^.^^" . (1.7) 

The vector field u is hypersurface forming if a; = 0. 

It is useful to define a representative length I along the worldlines of u, describing the volume expansion 
(contraction) behavior of the congruence completely, by the equation 



i0. (1.8) 
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^At this point there exists a lot of confusion in the literature. We think that in general it would be sensible to follow 
the conventions of Misner, Thorne and Wheeler [3], such that 770123 = \/~9 > r]°^^^ = — l/y'—g ; hence T)^vpa is the natural 
canonical 4-form representing the future pointing, righthanded unit volume element in four dimensions. Then — r]ixvp<T u" ~ 

T)apvp translates into the unit 3- volume element tyfe, and the decomposition of the gradient of the 4- velocity u reads 
^ pUjj — — Up-iiy + api, + ^ hpi, + ujpi, . Overall a minus sign would occur in all subsequent definitions, which involve rjpupa 
or its orthonormal frame form rjabcd '^•^ ■ Furthermore, applying these conventions, in all equations terms linear in the vorticity 
would come with the opposite sign. 
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1 . (1.9) 



It is also useful to define the Hubble parameter H and the dimensionless (cosmological) deceleration param- 
eter q: 

^-3°' ^- Kv,£)2 [e 

Work on the covariant 1+3 splitting of fluid spacetime geometries [1] was initiated first by Eisenhart and 
Synge and continued by Godel, Raychaudhuri, and the group comprized of Schiicking, Ehlers, Sachs and 
Triimper (see e.g. Ref. [4]). Extensions are due to, besides others, Hawking [5]. Nice reviews of this 
approach with particular focus on applications in relativistic cosmology were subsequently given by Ellis 
[6, 7]. 

The outline of the paper is the following: In section 2 we present and discuss the basic equations in the 
1+3 ortlionormal frame approach. General features of this approach and various applications have been 
discussed in the past by, for example, Pirani [8], Ellis [9] and MacCallum [10, 11]. Other useful references 
addressing general properties of the orthonormal frame approach are the books by Wald [12] and de Felice 
and Clarke [13] and the paper by Edgar [14]. The equations are given explicitly in terms of irreducibly 
decomposed quantities, including the Bianchi identities, which had not been given fully expanded in the 
aforementioned references. As the curvature is associated with, in principle, physical observations, it is of 
advantage to consider curvature related quantities as variables, as we discuss in this section. The present 
formulation is completely analogous to the Newman-Penrose approach [15]; we only replace a null congruence 
by a timelike congruence. It is also shown how one can naturally introduce local coordinates adapted to the 
existence of a preferred timelike vector field in the context of the orthonormal frame approach, following the 
work by Jantzen et al [16]. In section 3 we give a dimensionless formulation of the 1+3 orthonormal frame 
equations based on "expansion-normalized" variables. This is a direct generalization of the formulations of 
Wainwriglit and coworkers in the context of Bianchi cosmology and models with two commuting spacelike 
Killing vector fields (see e.g., Ref. [2] ). In section 4 we give a number of examples showing the usefulness 
of the formahsm developed in sections 2 and 3. First we show how one can include magnetic Maxwell fields 
in the so-called "silent" cosmological models [17] - [19]. We also discuss further possible generalizations of 
this class and the connection between the "silent" models and the algebraically special spacetimes in the 
context of the Newman-Penrose formalism. Thereafter it is shown how one can easily obtain a useful form 
for the equations of locally rotationally symmetric models (LRS) [9, 20] and spatially homogeneous models 
[21, 22], simply by deleting terms from the equations given in sections 2 and 3. We finally conclude with 
some remarks in section 5. 

We use the following index conventions for tensors: covariant spacetime indices are denoted by letters 
from the second half of the greek alphabet ( /x, z/, p, . . . = — 3 ), with spatial coordinate indices symbolized 
by letters from the second half of the latin alphabet ( i, j, fc, . . . = 1 — 3 ); orthonormal frame spacetime 
indices are denoted by letters from the first half of the latin alphabet ( a, fo, c, . . . = — 3 ), with spatial 
frame indices chosen from the first half of the greek alphabet ( a, /3, 7, . . . = 1 — 3 ). 



2 The 1+3 Orthonormal Frame Approach 

In the orthonormal frame approach one chooses at each point of the spacetime manifold a set of four linearly 
independent 1-forms { } such that the line element can locally be expressed as 

ds'' =riab<^''<^^ , (2.1) 

where rjab = diag [—1,1,1,1], ( = 1 ) , is a constant Minkowskian frame metric. The vectors { } 
dual to the 1-forms { w° } satisfy the relation 

{u>^e^,) = <5%. (2.2) 

In the 1+3 orthonormal frame approach to fluid spacetime geometries it is customary to choose the timelike 
frame vector Cq to be the 4-velocity of the matter fluid flow, although other choices are possible as well and 
will be mentioned later on. Along these lines we are now going to make a 1 + 3 split of the commutator 
relations as well as of the curvature variables and their field equations, part of which are constituted by the 
Bianchi identities. 
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2.1 The commutators 

The commutation functions, ^%c, are defined by 

[eQ,e;,] = 7^0;, Be , (2.3) 

where the frame vectors are understood to act as differential operators, ea{T), on any geometrical objects 
T. The covariant derivative is given by 

V.y^ = e„(y'^) + n„ , VbVa = e„(K) - T'^ab K . (2.4) 

If one assumes that there is no torsion and that the connection is related to the metric through 

VcVab = = - r'*,, r]db - r]ad , (2.5) 

one can express the Ricci rotation coefficients, r%c, in terms of the commutation functions by the relation 

Tabc = U^ad l\b + r]bd l\c - Vcd l\a ] ^ 7%c = " [ ^^c " ^^b ] , (2.6) 

where Tabc = iladT\^ and Ti^ai))c = 0. The latter property arises as a consequence of Eq. (2.5). 

Aligning the timelike direction of the orthonormal frame with the tangent of the preferred timelike 
congruence, sq = u {u°' = 6%, Ua = — S^a), the commutation functions with one or two indices equal to zero 
can be expressed in terms of the frame components of the kinematic quantities associated with the timehke 
congruence as defined in Eq. (1.4), and the quantity 

-In'^i^^et-e^Ud , (2.7) 

(where Gq := u''V ), which can be interpreted as the local angular velocity of the (to be chosen) spatial 
frame { Gq, } with respect to a second spatial frame { Sq, }, which is Fermi-propagated along gq = u. Secondly, 
the purely spatial components, 7^/37, are decomposed, following work by Schiicking, Kundt and Behr (see 
e.g. Ref. [21] ), into an object Oq, and a symmetric object riap as follows: 

-I"0^:=2a^p6%^+e0^,n'" . (2.8) 

eap^ is the totally antisymmetric 3-D permutation tensor with £123 = 1 = e^^^. The commutators are 
described by the expressions for ^%c- Their 1+3 decomposition leads to 

[eo,e„] = «„Go- [ |0(5^„ + (T^„ + e^„7(w^-fi') ] G^ , (2.9) 

[ea,G^] = - 2 ea/37 w'^ Go + [ 2 a[„ (5'^^] + eQ,/35 n'''^ ] g^ . (2.10) 

As an aside we here mention the conditions that a particular spatial frame vector be hypersurface orthog- 
onal (HSO) [12]. They are given by 

= (JaP + taPi ( - fiT ) O ^ 8 ^ 7 , 

= riaa no summation . 

2.2 The curvature 

The relation between the Riemann curvature tensor and the Ricci rotation coefficients is given by [2] 

R^cd '■= ec(r"6d) — BdiT^c) + T%d — T^d T%c — T%e 1%d , (2-12) 

while the Ricci curvature tensor and Ricci curvature scalar are defined by 

Rbd '■= R\ad , R '■= R'^^ab ■ (2-13) 
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The definition of the Riemann curvature tensor and the symmetry property T(^ab)c = imply Rabcd = — Rabdc 
and Rabcd = — Rbacd- The 16 Jacobi identities 

[[ea,ei],e,] + [[ei,e,],ea] + [[ec,ea],ei] = <^ e^^ij^c]) + [ab l'^ c]e = ^ , (2-14) 

correspond to the 16 equations 

R'^lbcd] = , Rabcd = Rcdab , (2-15) 

where the anticyclic relation R"'[i)cd\ = 0, together with the antisymmetric properties of the first and last 
index pairs of the curvature tensor, leads to Rabcd = Rcdab- 

The Einstein field equations with a nonzero cosmological constant can be written as 

i?% = T%-iT<5% + A<5% . (2.16) 

Here (and in the remainder of this paper) we use geometrized units, c = 1 = SttG/c^ [3], such that the 
dimensions of all dynamical variables can be expressed as integer powers of the single remaining nontrivial 
dimension [length]. 

In the phenomenological fiuid description of a general matter source the standard decomposition of the 
stress-energy-momentum tensor Tab with respect to a timelike vector field u is given by 

Tab = HUaUi, + 2q(aUb) +phab+T^ab ■ (2.17) 

Here the following matter fields arise: /i denotes the total energy density scalar, p the isotropic pressure 
scalar, q°' the energy current density vector, and TTab the anisotropic pressure tensor. We have that 

QaU" = , TTabu'' = , TT^a = , TTab = T^ba , (2.18) 

and the matter fields need to be linked by an appropriate thermodynamic equation of state in order to 
provide a coherent picture of the physics underlying a fiuid spacetime geometry scenario. 
The completely tracefree Weyl conformal curvature tensor is defined by [2] 

C^'cd := R^'cd - 2 6\ R%] + i n S% S",] . (2.19) 

When there exists a preferred timelike vector field u, it is also useful to decompose the Weyl tensor into its 
"electric part" 

Eac := Cabcdu'' u\ (2.20) 

and its "magnetic part" 

Hac--=*Cabcdu''u'' , (2.21) 

where the dual is defined by 

*Cabcd '■= ^Vab^^ defcd ■ (2.22) 

The "electric" and "magnetic" tensors are symmetric and tracefree and satisfy EabU^ = = Habu''. It 
follows that Eab = = Hat <=> Cabcd = 0. These definitions lead to [6] 

C^'cd = [ 4 S% 5% - r^'^Kf n'\d ] ug i?4 - 2 [ n^'ef S\c + 8K 8% n'\d ] ug H''h 

(2.23) 

Then, from Eq. (2.19), using Eqs. (2.16), (2.17) and (2.23), we obtain for the Riemann curvature tensor 
the expression 

R^\d = [ 4 6K 5% 6\ - n^\f n'\d ] ug Efn - 2 [ r,-\f 5% 5\r^ + 5^ n'\d ] Ug Rf^ 

+ 2 6% [{ix + p) «d] + g"! + g^j + ] + | (/^ + A) 6\ 6\ . (2.24) 

Inserting expression (2.12) for R'^^'cd into the lefthand side of Eq. (2.24) implies that one can obtain the 
Einstein field equations (by contraction on the indices a and c), the 16 Jacobi identities and also expressions 
for Eab and Hab in terms of the frame derivatives Bq and the basic variables 

{@,Ua,(^aP,(^a,^a,aa,nap,fJ',P,qa,''^ap} ■ (2.25) 
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The field equations 

eo(e) = -i02 + (e„+w„-2a„)(w")-2cr2+2a;2-i(/i + 3p)+A , (2.26) 

- i (5"^ [ (e^ +u^ + a^) (mT) + 2uj^W] + e^^^" [ 2 Q.^ cr^^^ - n^'^ us ] , (2.27) 
= ^l-\Q'^ +0"^ -J^ -2wa^" -\*R + K , (2.28) 

= g"-|5"^e^(0) + (e^-3a3)(a«'')+n"3^^^ 

- e"^^ [ (e^ + 2 - a^) (a;^) + n^s A ] > (2-29) 

*Sai3 ■■= e(„(a^)) + 6a/3 - i (5„/3 [ e^(a'^) + 6\ ] - (e|^| - 2a|^|) (n^)^) , (2.30) 

*i? := 2(2e„-3a„)(a")-ifo"„ , (2.31) 



where 



2na^n"'p — ri^^nap ■ (2.32) 



Equation (2.26), which is the (OO)-part of Eq. (2.16), is commonly called the Raychaudhuri equation. The 
symmetric tracefree (Q;/3)-part of Eq. (2.16) corresponds to the shear evolution Eq. (2.27), while Eq. (2.28) 
constitutes a generahzed Friedmann equation. Finally, Eq. (2.29) is the (OQ;)-part of Eq. (2.16). 

The objects *Sai3 and *R are the tracefree part and the trace of *Rai3, which is just the intrinsic 3-Ricci 
curvature of spacehke 3-surfaces orthogonal to the matter flow tangent u, if = [10]. However, when 

^ 0, then *Rap is not even a tensor; it is just a symbol defined by its relation to the objects Cq, and 11^/3. 
There are a number of ways one can define a "3-curvature" on the local rest 3-spaces orthogonal to u for 
the case when lo°' ^ 0, but we will not discuss these here. Instead we refer to Ref. [23]. 

The Jacobi identities 

eo(a") = -\ {8"^ ep + u" + a") (6) + ^ {ep + up - 2 ap) {a"^) 

- i e"^^ [ep +up-2 ap) (t^^ - Q^) , (2.33) 

eo(n"^) = - i G n"^ - {S^^" + 4") (w^' - fJ^') + 2 a^"^ n^''^ + (e^ + w^) (w^ - W) 

- [ (e^ + u^) (a'^^s) - 2 n''\ {u;s - ^a) ] , (2.34) 

eo(a;") = - ^ & uj" + a'^p cof^ + ^ n'^p u'^ - e'^'^'' [ ^ {ep - ap) {ii^) + cup fl^ ] , (2.35) 

= {ep-2ap){n"'^)-lQuj" -2(7"puj'^ + e"'^'' [ep{a^) + 2ujpn^] , (2.36) 

= (e„-«„-2a„)(a;") . (2.37) 
The "electric" and "magnetic parts" of the Weyl tensor 

Converting the covariant equation (4.16) in Ref. [6] into a relation with respect to an orthonormal frame 
gives 

- ^Sap [ (e^ +u.y+ ay) {u?) -20"^ -uP]+ e'\a [ 2 ^\y\ (Tp)s - np)^ us ] . (2.38) 
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Then, after combining this relation with the field equation (2.27), one easily derives 

EaP + ^ TfaP = \Q (^af} — (^ai CT^P - W/J - 2 ^(q, f2^) + ^ 6ap [ 2 Cj'^ + u'^ + 2 LOj fl'^ ] + *Sap . 

(2.39) 

From this expression it can be seen that the "electric part" of the Weyl tensor is closely related to *Sap- 

To obtain an explicitly tracefree equation for the "magnetic part" of the Weyl tensor one can combine the 
covariant equation (4.19) in Ref. [6] with the Jacobi identity ^w" = ii^uj^. This leads to the covariant 
expression 

H^i, = 2 oj,,) + h''(f, h^^-iVpOJa - I [2upUj'' + W^V pu" ] 

+ h%h\)npr.x{V^a^)u^ , (2.40) 

which, when formulated with respect to an orthonormal frame, gives 

Hap = (e(„ + 2m(„ + a(„) (w^)) + |nXyCT„/3 - 3n''(„o-^)^ - i^a^ [ {e^ + 2u^ + a^) {u}'^) - Sn^s (^'^^ ] 
+ a\',\){(^P)6) -iip)^u}s] ■ (2.41) 

Concluding this subsection, we briefiy review the deviation equation for timelike geodesies. This equation 
gives a direct measure of the focusing or diverging of two initially parallel fiducial timelike geodesies due to 
gravitational effects, as one follows along one of them towards say increasing values of its affine parameter. 
It was discussed in some length by, for example, Szekeres [24] ( see also the earlier discussion given by Synge 
and Schild [25] ) . The separation of the geodesies is encoded in a spacoliko vector field ^ ( using frame 
components: > 0), orthogonal to the tangent u (^qM° =0) of the reference geodesic (with ii" = 0). 

In principle the latter represents a potential observer. Then, on using the decomposition of the Riemann 
curvature tensor given in Eq. (2.24), the acceleration of the spacehke vector field ^ is expressed by 

= - [£%(u)-i7r%(u)] ^''-i[Mu)+3p(u)] r + Hr . (2.42) 

Supposing that ^ represents three eigendirections of a reference volume forming body, this relation provides 
a nice illustration of the different possible distortions the body may experience as it moves geodesically in the 
presence of nonzero spacetime curvature; specifically as measured in the local rest 3-spaces orthogonal to u. 
In these rest spaces three different effects contributing linearly to the acceleration of ^ can be distinguished: 
(i) the volume-preserving, shearing effects characteristic of tidal forces as exerted by the "electric part" of the 
Weyl tensor. Eat, and the anisotropic pressure of the matter fiuid, TTab (Note that TTat acts in diametrically 
opposed directions to Eat and thus tries to balance the shearing action of the latter.), (ii) the volume- 
shrinking effects of the total energy density, fi, and scalar pressure, p, (given that (// + 3p) > 0), and (iii) 
the volume-infiating, antigravitating effects induced by a nonzero, positive cosmological constant, A (which 
is an invariant spacetime constant and thus does not depend on the choice of u ) . The remaining information 
on the components of the Riemann curvature tensor as expressed by Eq. (2.24), namely the values of iJaf, 
and , can be found from measurements with respect to two further orthonormal frames { So } and { } 
that move in different directions relative to { }, which is ahgned with the tangent u [24, 3]. 

2.3 The Bianchi identities 

The Bianchi identities are given by [26] 

V[QJ?;,c]de = '^dCab''' = -'^ [aR%\ - \ ^\a'^ b\R = -'^ [aT%\ - \ 6\a'^ b\T ■ (2.43) 

We are now going to express these relations in terms of Eap, Hap and the irreducible parts of Tab- 
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Bianchi identities for the Weyl tensor 

Using the decomposition of the W^yl curvature tensor in terms of its "electric" and "magnetic parts" in the 
above expression (2.43) and subtracting the twice contracted Bianchi identities (2.48) and (2.49) (see below) 
where appropriate, leads to 

eo(S"^ + ^tt"^) = -i(/i+p)(T"^-e(B"^ + i7r"^)-i((5T("e^ + 2M(" + a(")(g^)) 

+ 3 cr("^ {E'^^'< - i 7r^)T) + i n\ i?"^ - 3 n^"^ i?^''^ 
+ [^{e^ + 2u^ + a^){q'^)-3a^s{E"'^ -^7r'^^) + 3n^sH"'^] 

- (u;^-2 n^) (E^\ + i 7r^\) + i n% qs ] , (2.44) 

eoiH"'^) = -QH"^ + ?,a%H^^^ -Ito^'^q^^ -\n^^{E"'^ -\t:"'^) 
+ 3n("^(£;'')^-i7r^5"') 

- \ qs + {to^-2 n^) ] , (2.45) 

= {ep-3ai3){E''^ + ^^7r^^)-lS"^ep{^i) + lQq"-^a^pq^ + 3u;pH"^ 

- e"^^ [ aps H% + lLOpq^+ nps (E^ + i Tr-*^) ] , (2.46) 

= (e^-3a^)(/?«^)-(/i + p)a;«-3t^^(£«^-i7r«^)-in«^/ 

+ [ i (e/3 - ap) (q^) + aps {E% + \ 7:%) - np, H% ] . (2.47) 

Bianchi identities for the source terms 

These arise after double contraction of Eq. (2.43), projecting the resultant relation respectively parallel and 
orthogonal to u. We get 

eo(/i) = - (/i + p)0-(e„ + 2M„-2a„)(g«)-(7„^7r«^ , (2.48) 

eo(9") = -l@q" -S"/'epip)-ifi+p)u'=' -{ep+up-3ap)iw"^)-a"pq'' 

+ e"^^ [ (ujp + np) q^ + nps tt% ] . (2.49) 



If the fluid flow congruence is irrotational, lo^ = 0, it is convenient to define a tensor encoding the conformal 
curvature properties of the 3-surfaces orthogonal to u. Throughout the literature it is mainly the conformally 
invariant symmetric tracefree 3-Cotton-York tensor (density) [27, 28] which is employed for these purposes. 
This tensor is divergence-free. We present it defined as the spatial rotation of the symmetric tracefree 3-Ricci 
tensor, a form due to Jantzen et al [23], which is equivalent to its usual definition. The symmetric tracefree 
3-Ricci tensor can be obtained by use of the Gaufi embedding equation [29]. In covariant form one gets: 

*S^v = {E,,^ + \ TTfj^) - I (Tp^ + a,jp a''ij - I (T^ h,j^ . (2.50) 

Then it follows for the covariant components of the 3-Cotton-York tensor ( h denotes the determinant of 

*C^, := /i^/'/i%/i^)/?p.«a(V"5^)m^ 
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= - h}'"" [ h% u^VrHp, + I If^, - 3 a% H,)p + o", H% h^, ] 

- 2u^ E\ + \(T\q'' ] , (2.51) 

where we have made a substitution in terms of the covariant form of Eq. (2.45) [6]. Again, specializing this 
expression to an orthonormal frame we arive at 

*Cap = e'^V (e|^| - a|^|) i*Sp)s) - Sn\a *Sp)j + \ n\ *Sai3 + Sap n^s *S^^ (2.52) 
= - eo{Hap) - I Hap + 3 (t'^(„ -ff/3)7 - 3 n^(a (^jS + \ rf-i cFaS cr^p + 2 Uap 

- 3 7r/3)^ + \ rC^ iTap 

— Sap [ (Tjs H''^ — n^s cr^e cr% + — n-^s tt'*'^ ] 

+ [ (e|7| - a|7|) (cr'/3) (T^e + 7r/3)5) - i e|^|(0) (T^)a - 2m|^I 

+ ^ap)^qs + 2n\^\Hp^s ] , (2.53) 
where Eq. (2.30) for *Sap is understood to be inserted into Eq. (2.52). 

2.4 The source 

There exists an abundance of possible choices for the matter source of the gravitational field. Here we will, 
however, consider only perfect fluids and Maxwell vacuum flelds. 

2.4.1 Perfect fluids 

A perfect fluid is described by the stress-energy-momentum tensor 

Tab = flUaUb + Phab , (2.54) 

where u°' is the fluid 4-velocity, fi the total energy density, p the pressure and h°'b the orthogonal projection 
tensor associated with u". It thus has vanishing energy current density (g" = 0) and anisotropic pressure 
{■JTab = 0). If we choose to take the comoving approach (eo = u), the twice contracted Bianchi identities 
(2.48) and (2.49) reduce to 

eo(/i) = -(/i+p)0, (2.55) 
= e^ip) + {p. + p) iia ■ (2.56) 

However, one does not always want to choose eo to coincide with the fluid 4-vclocity. In the case of several 
perfect fluids with different u one can only adapt to one of them. One might also have some other structure 
one wants to adapt to (e.g., in the spatially homogeneous models below we will choose eo to be the unit 
normal of the spatially homogeneous 3-surfaces). If the fluid is rotating one might want to use a 3 -I- 1 
initial value formulation and choose eo to be the unit normal associated with a family of spacelike 3- 
slices. We say that the fluid is tilted if u ^ Bq. In this case we obtain a stress-energy- momentum tensor 

Tab = I^UaUb + 2 q(a Ub) + P hab + T^ab with 



2 (r. .. _ .. 

(2.57) 



/i = V^{fi + p)-p, p = j{p+p)T^v^ +p , 

q'' = r2(/i-|-p)w% TTab = r'^{fj+p){VaVb-lv'^hab) , 



where = F {u'^ + v'^), Ua v"" = {),V := l/Vl — v"^, '■= Va v"'. For a linear barotropic equation of state of 
the form p{ji) = (7 — 1) /« these equations reduce to 

/i = r^Gfi, p = iG-i[^3-27)^;2+3(7-l)]/i, .3 58) 

= 'yG ^V^fl, TTab = -lG ^ {VaVb- \v'^ hab) , 

where G := 1 4- (7 - 1) w^. 
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The source equations (2.48) and (2.49) yield evolution equations for /x and v"' respectively. The one for 
/i is obtained directly by inserting the above expressions into the /x-equation for a general stress-energy- 
momentum tensor. The one for v"' is quite complicated and we will only show how one can obtain the 
Pifi) = (7 ~ 1) A* case by taking the following linear combination of the /i- and g^-equations: 

eo(^") = [ 1 _ _ 1) ^2 ] [ - 7^" eo(/i) + [ 1 - (7 - 1) ] eo(g") + 2 (7 - 1) ^" vp eo(g^) ] . (2.59) 

This equation will be evaluated exphcitly in dimensionless form for spatially homogeneous perfect fluid 
models below (it is quite cumbersome in the general spatially inhomogeneous case). One can obtain an 
evolution equation for // once one has one for v". However, it is often preferable to keep /x and use the 
constraint equation (2.28) to solve for /x in terms of other variables. 

2.4.2 Maxwell vacuum fields 

An electromagnetic field is characterized by the antisymmetric field strength tensor [7] 

Fab=UaEb-UbEa+riabcdH'u'^ , (2.60) 

where Ea denotes the electric field and Ha the magnetic field respectively and have the property = 
Eau" = Hau". Treating the Maxwell vacuum as a fluid with 4- velocity u the contributions to its stress- 
energy-momentum tensor in an orthonormal frame can be expressed by 

fi = ^{E^E" +H„H") = 3p, (2.61) 
q" = e"^'> Eg H., , (2.62) 
7r„^ = -E^Ep-H^Hp + \5^p{E^E^ + H^H^) . (2.63) 

The sourcefree Maxwell equations with respect to this orthonormal frame are given by 

eo(£") = -lQE" + a"pE^ -n"pH^ +6"^-^ [{ep+up-ap){H^)-{ujp-np)E^] , (2.64) 

eo(Jf") = -l&H'' + a''pH'^ + n''pE'^ -e'"^^ [ie(3+up-ap)iE^) + iu;p-n(3)H^] , (2.65) 

= (e„-2a„)(E«) + 2a;„Jf" , (2.66) 



= (e„-2a„)(iJ")-2tj„£" . (2.67) 
The covariant form of these equations was discussed in Ref . [7] . 

2.5 Introducing local coordinates 

For many purposes one does not need to introduce local coordinates; it is often not even advisable to do 
so. However, many things do rely on the introduction of a local coordinate system, and if one wants to 
introduce coordinates one might as well do so in a geometric way reflecting the problem at hand. In the 
present context we have a preferred timelike congruence u. The natural way of introducing coordinates 
adapted to this particular structure is the 1+3 threading approach, recently discussed by Jantzen et al [16] 
and Boersma and Dray [30]. In this approach one expresses the orthonormal frame as follows: 

eo = dt, e„ = ej (Mi dt + 5,) . (2.68) 

This leads to the following expression for the 4-geometry 

^^'g^^ = - 60 60 + ^"'^ e„ eg) eg 

= - dt®dt + (5"^ e„' ep^ {Mi dt + dt) {Mj dt + dj) , (2.69) 
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where M = M(t,x^) is the threading lapse function and Midx^ = Mi(t,x^) dx^ is the threading shift 1-form 
[16, 30]. The triad vector components e„* = e„*(f, a;*) describe the threading metric 7'-' = 6"^ ea ep^ ( this is 
just the spatial projection tensor expressed in local coordinates; for a discussion on its physical significance, 
see Refs. [16], [30] and [31] ). Thus, in the present 1+3 threading approach, i is a coordinate along the 
timelike congruence, while the parametrize the difTerent flow lines of the timelike congruence. This is 
the "dual" formulation of the more familiar 3+1 decomposition. In that formulation one imposes a causal 
condition on 3-surfaces, which are assumed to be spacelike, while no condition is imposed on the flow lines 
threading the 3-surfaces. In the threading approach one imposes a causal condition on the congruence, which 
is assumed to be timelike, while no causal condition is imposed on the 3-surfaces. 

Note that M,; ^ if one has a rotating congruence, uj" ^ 0, which is the case when the local rest 3-spaces 
are not hypersurface forming. Of course, one still has a considerable freedom in choosing the threading 
shift Mi; a freedom which just corresponds to the freedom of choosing a foliation such that the timelike 
congruence is nowhere tangent to it. 

Expressing the commutators (2.9) and (2.10) in terms of the variables M, Mi and yields the relations 



ej [ {Mi),t + {M,i + Mi M,t) ] = m„ , (2.70) 

{ej),t = -Me^' [ i0<5^„ + a^„+e^„^(c.'^-n'^) ] , (2.71) 

M e[J epf (Mi J + Mi,t Mj) = e^p^ u;^ , (2.72) 

2 e^J [ (e^]^),i + (e^]^),t Mi ] e^^ = 2 a[„ <5^^] + e^ps , (2.73) 
where := dtf, f,i := dif, and where e"i is defined through the relation 

e"ieJ=Pi. (2.74) 



Note that one can insert the expression for (eQ,'),( given in Eq. (2.71) in Eq. (2.73), thus obtaining an 
expression involving only spatial derivatives. 

As an example of the threading approach we consider Taub's comoving coordinates for isentropic perfect 
fiuids with equation of state p = p{p) [32]. Taub's comoving coordinates correspond to the fohowing threading 
lapse function and shift 1-form: 

M = r-\ r := r, exp r — , Mi = Mi(x^ ,x'^ ,x'^) , = r S^o , (2.75) 

where r is related to the fiuid's enthalpy density {fl+ p). With this choice of threading lapse function and 
shift 1-form the momentum conservation equation (2.56) will be automatically satisfied, as can also be seen 
from Eq. (2.70). 

If the congruence is nonrotating (w" = 0) , it naturahy gives rise to a fohation of the spacetime manifold 
by spacelike 3-shces. In this case one can choose to focus on the slices rather than the congruence. This 
situation suggests to take the 3 + 1 slicing point of view in which one introduces local spatial coordinates 
on the slices, while one introduces timelines that do not, in general, follow the original timelike congruence 
(see e.g. Refs. [16] and [30], as well as the standard texts in Misner, Thorne and Wheeler [3] and by York 
[33] ). Here we will choose to express the spatial triad { } in a noncoordinate basis, since this is useful 
in the context of, for example, Bianchi cosmology (a similar formulation could have been introduced in the 
threading approach; but we will not pursue this any further here). In the slicing approach the orthonormal 
frame is expressed as follows (Note that, deviating from our conventions, here i, j, k are also used as spatial 
indices for the noncoordinate basis.): 

eo = (dt -N'ei) , e„ = ej , (2.76) 
which leads to the following 4-geomctry 

(4'g-i = - eo « eo + (5"^ e„ (X) 

= - N-'^ (dt - N' Bi) (dt - Sj) + (5"^ eJ ep^ , (2.77) 
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where TV = N{t,x'^) is the lapse function and TV* = N^(t,x^) is the shift vector. The inverse e"i{t,x'') to 
ea-'{t,x'^) (e"j Bq,-^ = S^i) yields the spatial metric tensor g,j = Sape"ie^j on the spacelike 3-slices expressed 
in the basis = ei^{x'') dj. 

The above basis yields the commutation relations 

[dt,ei] = , [ei,ej]=^%{x')ek , (2.78) 

which, when compared to the commutators (2.9) and (2.10) and setting = 0, leads to 

iV-i eJei{N) = M„ , (2.79) 

- Nep' [led^a+ - e^„7 ] , (2-80) 

e\ [ 2e[Jei{ep]'') + e[Jepfj% ] = ^a^^S^'pj + e^psn''' . (2.81) 

Geodesic and nonrotating congruences ( = = ) imply that one can introduce a synchronous local co- 
ordinate system ( TV = 1, TV* = 0; M = 1, = ). These congruences have several interesting consequences. 
Equations (1.9) and (2.26) imply 

-l = ±[2a' + ^{^l + 3p)-A] . (2.82) 

Thus (/i + 3p) > and A < imply that q > 0. and that a singularity forms in the past (future) with 1 = 
if > (0 < 0). The singularity may be a curvature singularity or just a crushing singularity (see e.g. de 
Felice and Clarke [13] and MacCallum [10] ). Equation (2.28) and the above expression for q lead to 

9 = 2-^ [ i(/^-p)+3A-*i?] . (2.83) 

Thus if A > 0, /i > p, *J? < 0, this leads to the inequality q <2. Hence, if A = 0, (/i + 3p) > 0, (n - p) > 0, 
*R < 0, it follows that < <? < 2. 

Of course, there are other ways local coordinates can be introduced. For example, one might want to 
focus on events in the neighborhood of a single line of the timelike congruence one is interested in. In this 
case it is natural to introduce Fermi normal coordinates (for a discussion see e.g. Ref. [3] ). 

2.6 Evolution equations and constraints 

The usual initial value problem in General Relativity is associated with a 3 + 1 slicing formulation. The 1+3 
approach, with its associated threading formulation, is uncharted territory ( for a discussion see [23] ). We will 
thus instead discuss the 3 + 1 approach in the present tetrad formulation. In the 3 + 1 initial value formulation 
one usually uses the spatial metric tensor of the 3-slices and its extrinsic curvature tensor as variables. The 
lapse function represents the freedom in choosing a spatial foliation (choice of time) , while the shift vector 
represents the freedom in choosing a threading (choice of spatial gauge) . In the present formulation the triad 
components, Bq,*, replace the spatial metric {gtj = Sape^ie^j => 5ape'^[ie^j] = 0), while the expansion 
scalar and the shear tensor correspond to the trace and tracefree part of the extrinsic curvature tensor 
respectively; the frame derivative gq is just the normal derivative to the spacelike 3-sHces. The role of 
derivatives in the tetrad formulation becomes trivial, since all tetrad components are spacetime scalars. 
Note that it is only gq that contains a temporal partial derivative (gq = N~'^ {dt — TV'cj) , = ea'si). 
Thus one obtains an evolution equation in standard form when one has a eo-derivativc on the lefthand side 
simply by moving the shift part to the righthand side. Thus equations involving eo will be referred to as 
evolution equations, while equations not involving eo are said to be constraint equations in analogy to the 
usage of EUis in the covariant formulation [6, 7] (even though the interpretation in the rotating 1+3 case 
is not clear). 
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It is seen that by extending the number of variables, from the triad, expansion and shear variables to 
all commutation functions and curvature variables, one obtains additional evolution equations as well as 
constraint equations. However, there are no evolution equations directly obtained for the acceleration, ii", 
and the Fermi- rotation, fl". 

This is not surprising, since one can choose the timelike congruence so that one can obtain any convenient 
value for the acceleration u". This freedom is also seen when choosing local coordinates, in particular, one 
can choose any lapse function (recall Eq. 2.79). However, if one has a source which one wants to adapt to, 
this might provide an evolution equation for the acceleration. For example, in the comoving perfect fluid 
case the pressure gradient determines the acceleration. If one has a barotropic equation of state, p = 
one can obtain an evolution equation by using the commutators in conjunction with the source equations 
(2.48) and (2.49). This leads to 

eo(«a) = ea(|§e) + (|£ - i ) e - [ + e^„^ (o;^ - 0^) ] up , (2.84) 

which for a linear barotropic equation of state of the form p(/x) = (7 — l)fJ reduces to 

eo(wa) = (7 - 1) e„(0) + i (37 - 4) M„ - [ + e^„^ (to'' - n^*) ] up . (2.85) 

Neither is it surprising that one docs not have an evolution equation for the Fcrmi-rotation, f!", since one 
can choose 0" to be anything one wants (e.g., one can choose a Fermi-transported spatial frame, fl" = 0). 
Thus f2" represents the freedom of choosing the spatial frame { } and hence plays a role analogous to the 
shift vector, which represent the freedom in choosing spatial local coordinates. The Fermi-frame is not the 
only interesting choice for ft". Other possible choices are those which diagonalize a^p (shear eigenframes) 
or Eap. Of considerable interest is also the "co-rotating" frame, defined by fl" = u". 

We also note that we only have an evolution equation for the combination {E"^ + ^tt"^), and no one for 
the pressure p. To obtain individual evolution equations for these quantities requires further assumptions 
about the source. In the case of perfect fluids with, for example, p = p(/i) as an equation of state, the 
pressure p and TVap arc determined by the energy density // and the velocity v" ( sec Eq. (2.57) ) which have 
evolutions equations. For Maxwell fields one obtains fx, p, iTap from Ea and Ha (see Eqs. (2.61) - (2.63)) 
which also have evolution equations. Suggestions for a fully relativistic thermodynamic treatment of general 
viscous fluid matter sources with nonzero q" and iTap were given, for example, by Israel and Stewart [34, 35]. 

If one expresses the frame { } explicitly in terms of the lapse function, shift vector and the triad 
variables and calculates the Einstein field equations, one obtains a second-order system of partial differential 
equations. The commutation functions arc defined through the commutators, and the Jacobi identities are 
identically satisfied. However, one can instead view the commutation functions as a set of new independent 
variables. In this case one obtains a first-order system of equations in all derivatives described by the 
commutators, the Einstein field equations and the Jacobi identities, which are now "elevated" to nontrivial 
equations. The Ricci identities provide a definition for the Riemann curvature tensor. Using this definition 
one finds that the Bianchi identities are identically satisfied. 

However, the curvature is a directly measurable object. This makes it desirable to consider the curvature 
components as variables as well. Thus we can continue and view the curvature components as independent 
variables. In this case the Bianchi identities are "elevated" to nontrivial field equations. Note that the system 
of equations for tetrad variables, commutation functions and curvature variables is highly redundant (we 
could have stopped at, e.g., the commutation function level with the Jacobi identities and the Einstein field 
equations). From a physical point of view, in the present formulation it would be desirable to start at the 
bottom of the derivative ladder with the curvature and the Bianchi identities. Can we go "upwards" towards 
lower derivatives and only use, say curvature and kinematical quantities? Do these form a complete and 
closed system? In general the answer is, unfortunately, no. In general we always need some tetrad variables, 
commutation functions and commutator equations. The redundancy, in the general case, is described by the 
Papapetrou identities [36, 37] (for a pedagogical review see Ref. [14]). Thus one can pick out interesting 
subsets of variables and equations in the general case. For practical reasons we are usually interested in special 
cases, for example, in special Petrov types. Unfortunately these special cases need further consideration and 
one might have to consider higher derivatives for checking consistency. The "silent" models discussed below 
are examples of this. 
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However, it should be pointed out that one can take the opposite route and go "downwards" towards 
higher derivatives, since it is possible to start at the curvature (or commutation function) level and take 
a finite number of covariant derivatives and obtain a complete description of the geometry. This is the 
equivalence problem approach along the lines of work by Cartan [38] and Karlhede [39]. In this approach 
one obtains, step by step, the information contained in the tetrad and the commutators, and after a finite 
number of covariant derivatives one has obtained all the information needed to construct the full tetrad, that 
is, the metric [40]. 

2.7 Irreducible tracefree decomposition 

In the case of the covariantly defined spacclikc symmetric tracefree tensor fields Uaa, Eap, Hag and iVag, it is 
sometimes convenient to introduce a new set of variables adapted to the tracefree condition and the invariant 
quadratic forms , E"^ , H"^ and tt^. For example, for we define the irreducible frame components 

(T+:=-|(Tii = I ((722 + 0-33) , (T_ := ((T22 - 0-33) 

(2.86) 

fTi := -\/3(T23 , fT2 := -\/3(T3i , (T3 := -\/3(Ti2 , 

leading to 

cj^ = \[ {a+f + {a.f + {a^f + {a,f + {a,f ] , (2.87) 

where the definition is motivated by the desire to simplify the kinematic part of the /i-equation (2.28) 
as far as possible (this equation plays a key role in determining the dynamics of, for example, spatially 
homogeneous models). If one decomposes the shear tensor one should also similarly decompose the other 
spacelike traceless symmetric tensors as well. In principle one could extend this procedure to also include 
the spatial commutation functions tIq,^, that is, splitting them into a trace and a tracefree part, where the 
latter is further subdivided according to the scheme above. However, we will not take this step within this 
work as Uap usually does not have a covariant meaning in the 1+3 picture, but very much depends on the 
choice of the spatial frame { }. 

3 Dimensionless Formulation 

Dimensionless variables have at least a three decade long history in General Relativity. Bondi was perhaps 
the first to use scale-invariant variables in order to obtain a dimensional reduction of a problem. He did 
this for static spherically symmetric perfect fiuid models with p(/i) = (7 — 1) in 1964 [41] ( see also Ref. 
[42] ) . Independently Collins used the scale invariance of the Einstein field equations to produce reduced 
sets of equations in Bianchi cosmology [43] . Wainwright and collaborators have introduced the concept of 
"expansion-normalized" dimensionless variables in Bianchi cosmology and inhomogeneous models with two 
commuting spacelike Killing vector fields [44, 2] ( closely related variables have also been used by Rosquist 
and Jantzen [45] ) . Expansion-normalized dimensionless variables lead to relatively simple equations and are 
natural in a cosmological or gravitational collapse scenario ( for a recent review see Ref. [2] ) . On the other 
hand, in other situations, like for example for static models or for maximal slicings, one is forced to consider 
other possibilities. For special cases there might also be mathematical reasons for considering other choices 
(see e.g. Ref. [46]). However, for problems in cosmology and gravitational collapse, expansion-normalized 
dimensionless variables seem to be a natural first candidate, and this is how we will proceed. Here we 
will generalize the work by Hewitt and Wainwright [44] on comoving perfect fluids in a geometry with two 
spacelike commuting Killing vector fields to general timelike congruences and geometries. 
We first introduce expansion-normalized dimensionless diflFerential operators 

8. := ^ , (3.1) 
then expansion-normalized dimensionless commutation functions 

I Ua, S«/3, IFq, TVq,^ } '~ ''^ '^"'^ } / ® ' ('^•2) 
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and finally expansion-normalized dimensionless source and curvature variables 

{ fl, P, Qa,IiaP, ^A, £aP,'HaP } := 3 { H,p, ^a, 7r„/3, A, Eap, Haf} } / 0^ . (3.3) 

For the traceless decomposition variables of the shear tensor, Eq. (2.86), one can use the following definition, 
which is aimed at simplifying the 0-equation arising from Eq. (2.28): 

{S±,Ei}:={a±,a, }/©, (3.4) 

where j = 1, 2, 3, and from which it follows that 

y? := f = (S+)2 + (S_)2 + (Ei)2 + (S2)2 + (Sa)^ . (3.5) 

Then one can make a similar transformation for the other spacelike symmetric tracefree tensors: 

{ U±,ni,£±,£i,n±,ni } :=3 { 7r±,7ri,E±,Ei,H±,H, } / 9^ . (3.6) 

Since the only variable that carries dimension is 0, with dimension [ length ]^^, the equations associated 
with these variables will decouple from the remaining ones. The decoupled equations are given by 

Decoupled equations 

do® := -il + q)@ (3.7) 
a„0 := -r„0, (3.8) 

where q is the deceleration parameter defined in equation (1.9). These expansion-normalized dimensionless 
variables lead to the following expressions for the commutators, curvature equations, and Bianchi identities. 

3.1 The commutators 

[ao,a„] = - [r„-3f7„]ao+3 [ig<5^„-S^„-e^„^(iy^-iZ^)] , (3.9) 

[a„,a^] = -Ge^p^W^ do+[2ir[^+3A[^)6^p]+3e^psN^^] . (3.10) 

3.2 The curvature 
The field equations 



q 



{da-ra+SUa-6Aa)U" + 2'S'^ -GW"^ + Un + 3P) - Qa , (3.11) 



aoS"^ = (g-2)S"^-F(5T("a^-r("-l-3L/^("-H3A(")f/^5-H6M^("i?^)-Fn"^-5"^ 

+ 3 eT"'(" [ 2 Y,^\ - Us ] , (3.12) 

Q. = l-T,'^ + 3W'^ + 6WcR" -IC-flA , (3.13) 

= Q'^ + ^r" + {d0-r0-9A(3)i:"'^ +3N"pW'^ 

- e"^'' [(dp-rp + 6Up-3 Ap) + 3 Nps ^\ ] , (3.14) 
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where 

Saf} ■= 3 *Sa0 I 0^ 

= - r(J + 3 - i 5a0 [ {d^ - r^) A'>+3 B\ ] 

- e''%(a|^| -ri^i -6A|^|)A'>)a , (3.15) 

/C := - (3 / (2 e^) = - (2a„-2r„-9A„)A" + |S"„ , (3.16) 

Bap = 2 7V„^ N-'p - iV„^ , := W" . (3.17) 

The Jacobi identities 

doA" = gA" + ir"-L/" + i(a^-r^+3f73-6A^)S"^ 

- i e"^'' {dp -rp + 3Up-6Ap) {W^ - R^) , (3.18) 

SoiV"^ = giV"^ - (JT("a^ -r(" + 3J7(") (VF^^ -i?^^) + 6E("^A^^^^ 
+ 5"^ (a^ - + 3 J7^) {W - R^) 

- eT"'(" [ (a^ - + 3 Lf^) Y.'^\ - 6 NI^\ {W5 - Rs)\ , (3.19) 
doW" = {q-l)W'' + 3i:"pW'^ + ^N"pU'^ -e"'^'' {dp- rp- 3 Ap)U^ + 3WpR^^ , {3.20) 

= {dp-rp-6Ap)N"l^ -2W" -6'S"pW^ +e"'^^' [{dp-rp)A^+6WpR^] , (3.21) 

= (a„-r„-3f7„-6A„)W« . (3.22) 
The "electric" and "magnetic parts" of the Weyl tensor 

Sap = -(ao-g + l)S„^ + (a(„-r(„ + 3J7(„ + 3A(„)f7^)-3E„^S^^-3W^„VF^ + in„^ 
- ^Sc^p [ {d^ -r^ + 3U^ + 3A^)U'^ -GT,'^ -3W'^ ] 

+ 3 [ 2 J?|^| E^), - iV^)^ Us ] . (3.23) 
The field equation (3.12) combined with the above expression (3.23) for Sap leads to 

^aP + ^ ^aP = ^aP - 3 E - 3 Wa Wp — 6 W^^ Rp) 

+ \6^p[2T? + 3W'' + QW^R^]+Sa,p . (3.24) 

The dimensionless version of Eq. (2.41) is 

'Hap = (a(„-r(„ + 6f7(„+3A(„)W^) + |iV\S„^-9iV%S^)^ 

- I Sap [ (a^ - + 6 + 3 A^) -9 N^s S^'' ] 

+ e'\a [ (a,^! -ri^i -3A|^|)E/3)5-3iV^)^T^5 ] . (3.25) 
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3.3 The Bianchi identities 

Bianchi identities for the Weyl tensor 

do{S"^ + ^U"^) = (2g-l)f«^+ (g+i)n«^- f(n + p)E«^ 

- i - 2 r(" + 6 tl^" + 3 A^") Q^^ 

+ 9 S("^ (f - \ n^'T) + f 7^"^ - 9 iV("^ 

+ iJ"^ [i(a^-2r^+6?7^ + 3A^)Q^-9S^5(f^^-inT^)+9A^^5?^^''] 
+ eT"'(" [ (a^ - 2 + 6 f/^ - 3 A^) -H^^^ 

- 3 (W^^ - 2 R^) {Sl^\ + i n^)^) + I iV^)^ Qs ] , (3.26) 

aoi?"^ = (2g- 1)?^"^ + 9E("^?^^)t - f VF("Q^) - f A^\(f"^ - in"^) 
+ 9A^("^(f^'^-in^'^) 

- 3 J"^ [ E^a - i VF^ + iV^a (^^^ - ^ n^^) ] 

- e^^(" [ (a^ - 2 - 3 A^) {£^\ - \ + 6 

- l^^\Qs+2,{W^-2R^)'Hl^\] , (3.27) 



+ 9 VF^ - 3 e«^^ [ + lWfiQ^+ Nps {£% + ^ 11%) ] , (3.28) 

= {dp-2rp-9A0)'H''^ -3{n + P)W''-9Wp{£"^ -^U"^)-lN''pQ^ 

+ 3 e"^-^ [ i (a^, - 2 - 3 Ap) + ^ps {£% + i U%) - Hps n% ] . (3.29) 

Bianchi identities for the source terms 

doQ = i2q-l)n-SP-{da-2ra + QUa-QAa)Q" -ST.apU"'^ , (3.30) 

doQ" = 2iq-l)Q"-6"^idp-2rp)P-3in + P)U"-idp-2rp + 3Up-9Ap)U"^ 

- 3 S"^ + 3 e"^'>' [ {Wp + Rp) + Np^ U% ] . (3.31) 



Eqs. (2.53) and (2.52) convert into 
Cap '■= 9 *Cap I 0^ 

= e-^V (a,^! - 2 ri^i - 3 A|^|) Sp)s - 9 N\a Sp)^ + | iV^^ + 3 N^s S^' (3.32) 
= - (ao - 2 g + 2) -H^^ + 9 E^(„ ?^^)^ - 27 iV'^(„ S''^) E^^ + f TV^ E^^ ^^p + 6 iV„^ E^ 

- 9iV%n^)^ + f 7V\n„^ 

- 3 (5„/3 [ - 3 iV^5 ^^ S% + - 3 N^s H^"' ] 

+ 3e^% [ (a,^! -2r|^| -3A|^|)(S'^)S5, + in^),)+r|^|E^), -2J7|^|f^), 

+ i ^(3)7 '3^ + 2 -^ItI ^/3)^ ] • (3.33) 
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3.4 The source 

3.4.1 Perfect fluids 

For the comoving (cq = u) perfect fluid case we have: 

do^ = (2g-l)n-3P, (3.34) 
= ida-2r^)P + 3in + P)U^ . (3.35) 

The evolution equation for the acceleration , Eq. (2.85), now takes the form 

dot/a = - (7 - 1) ra + [ 3 (7 - 1) + g ] f^a - 3 [ + e^„^ (W'' - R^^) ] Up . (3.36) 

For a tilted perfect fluid we have 



n = r2 (n + p) - p , p = i (0 + p) r2 + p , 

Q'^ = r2(n + P)««, Uab = r^n + P)(VaVb-lv^hat) . 



(3.37) 



(3.38) 



For a linear barotropic equation of state of the form p{jl) = (7 — 1) /i these equations reduce to 
n = V^Cn, P = iG-M (3-27)^;2 + 3(7-l) ]n , 

where G := 1 + (7 - l)v'^. 

The source equations yield evolution equations for and v"". The one for fl is obtained directly by 
inserting the above expressions into the fJ-equation, Eq. (3.30). The one for v°' can be obtained for the 
p{jj) = (7 — 1) /i case by taking the following linear combination of the fJ- and Q'*-equations, Eqs. (3.30) 
and (3.31): 

dov" = ^ ^ ^ ^ _ 1^ _ J [ - 7^" + [ 1 - (7 - 1) ] doQ" + 2 (7 - 1) vp doQ^ 

+ 2il + q)[jnv" -il-ij-l)v^)Q" -2ij-l)v"vpQ^ ]] . (3.39) 
This equation will be evaluated explicitly for the spatially homogeneous perfect fluid models below. 

3.4.2 Maxv^rell vacuum fields 

We now introduce the following expansion-normalized dimensionlcss electric and magnetic field variables: 

{£a,na}:=V3{E^,Ha}/e . (3.40) 
This leads to the contributions to the Maxwell stress-energy-momentum tensor 

n = liEa^" +'Ha'H°') = 3P, (3.41) 

ga ^ t^'^^Epn^, (3.42) 

n„/3 = -£c£p-ncnp + lSc,p{£^£^ + n^n^) . (3.43) 

The sourcefree Maxwell equations are given by 

do£" = {q-l)£" + 3i:"p£'^ -3N"pn'^ 

+ [{dp-rp+SUp-S Ap) - 3 {Wp - Rp) £^ ] , (3.44) 

= iq-i)n" + 3i:"pn'^ + 3N"p£'^ 

- e"^-^ [{dp-rp+3Up-3Ap)£^ + 3{Wp-Rp)n^] , (3.45) 
= (a„-r„-6A„)f" + 6PF„?^" , (3.46) 
= (a„ -r„ -6A„)?^" -6W^„f" . (3.47) 
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3.5 Introducing local coordinates 

In the 1+3 threading approach the quantities M and Ca carry dimensions [length] and [length]"^ respec- 
tively, while Mi is dimensionless. The relation da = Ssq/Q leads to the introduction of the relative threading 
lapse function 

M:=MQ , (3.48) 
while the relation da = Se^/G leads to the introduction of the scaled triad 

Ej:=^, E"i:=Qe"i. (3.49) 

Thus we have 

eo = @M-^dt do = ^M-^dt, (3.50) 

e„ = e Ej (Mi dt + di) da = 3 Ea' {Mi dt + di) , (3.51) 

which leads to the possibility of expressing the 4-geometry as follows 

(4)g-i = [ - dt®dt + 6"^ Ej Ep^ (Mi dt + di) ® (Mj dt + dj) ] . (3.52) 

Expressing the "dimensionless commutators" in terms of A^, Mi and Bq,' leads to the dimensionless form of 
Eqs. (2.70) - (2.73): 

3Ej [{Mi)^t+M-^{M,i + MiM,t)] = -[r„-3J7„], (3.53) 

{Ej),t = MEp' [^S^'a-^^a-e^aAW -R'')] , (3.54) 
M E^' Epf {Mij + Mi,t Mj) = eap^ W , (3.55) 

2E[J [ {Epf),i + {Ep^^),tMi ] E-'j = 2 (A[„ + r[„)<5^^] + e^ps N'-* , (3.56) 

where one can insert the expression for {Ea^)^t given in Eq. (3.54) into Eq. (3.56), thus obtaining an 
expression only involving spatial derivatives. 

In the sHcing approach to the nonrotating {to" = 0) case we assign dimensions [length] and [ length ]^^ 

to the quantities N and e^', while the shift vector TV* and the spatial frame vectors are taken to be 
dimensionless ( the components of e.; are assumed to be given functions of the spatial coordinates , which 
themselves are taken to be dimensionless; for a discussion of coordinates and dimensions of spacetime tensors 
see Ref. [47] ). The relation do = 3eo/0 leads to the introduction of the relative lapse function 

JV:=NQ, (3.57) 

while the relation da = Ssa/Q gives 

Ej:='-^, E"i:=Qe"i. (3.58) 

Thus we have 

eo = QM-'^ (dt-N'ei) => do = SAf-'^ {dt -W Si) (3.59) 

e„ = QEa'ei da = 2,Ejei, (3.60) 

which leads to the 4-geometry 

(4)g-i = [ - 7^-2 (^Q^ _ jf ^ (^Q^ -N'ej) + S"^ Ej Ep^ ] . (3.61) 
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The relations (2.79) - (2.81) take the form: 



- [ra - SlJa] , 



(3.62) 



+ MEp' [^q8^a- S^a + e^a-, W ] , (3.63) 

-E'k[2E^Jei{Ep^^) + E^jEp^itij] = 2 (A[„ + r[„)<5-^^] + 7V^^ . (3.64) 

As mentioned earlier, one can introduce synchronous local coordinate systems {N = 1, = 0) for 
nonaccelerating and nonrotating congruences. With this choice of lapse the time variable tp measures proper 
time (also sometimes denoted as "clock time"). However, if one is so inclined, one can choose a new time 
variable characterized by TV = N{t). Note that individual curves of the congruence have constant spatial 
coordinates (since N' = 0). Thus one can choose N = ±3/0, A/" = ±3, on an individual curve (where 
+ is chosen when > and — is chosen when < 0). With this choice Eq. (3.59) leads to the relation 
do = ±c?r; where r is the time variable associated with the above lapse. From the definition of in terms 
of £ it fohows that d-j-i = ± 1 which in turn yields 

£ = £,e±(^-^-) . (3.65) 

Thus it follows that r — > =Foo at the singularity where £ — > 0.^ The equation SqQ = — {1 + q)Q yields 

@ = @^e'^0'+''^'^\ (3.66) 
for an individual curve of the congruence. From TV = ± 3/0 it follows that 

tpi*) - tpis) = Atp = ± r I dr = ± ^ r dT , (3.67) 

where tp{s) is the time at the singularity. This leads to 

i 0, Atp = H^Mp = ± j ' ■f-. ^^^"^ dT . (3.68) 

Hence the inequality < g < 2 leads to (dropping the *) 

i < HAtp < 1 , (3.69) 

where H Atp < 1 requires A < 0, ip +3p) > 0, while 1/3 < H At requires A > 0, {p - p) > 0, *R < 0, 
and thus both inequalities are satisfied if A = 0, (/x + 3p) > 0, {p — p) > 0, *R < along a given curve. 
The present discussion generalizes the results of Wainwright in [2] from orthogonal spatiahy homogeneous 
perfect fluid models to general nonaccelerating nonrotating congruences. As an example one can mention 
nonrotating inhomogeneous dust models. Current observations of the age of the Universe, Atp, and the 
present value of the Hubble parameter, H, suggest that HAtp may satisfy the inequality HAtp > 1 ( see e.g. 
Ref. [2] ). Thus, observations may cast doubt on these cosmologies in the near future. Apparently one may 
thus be forced to consider other sources, a positive cosmological constant, or rotation. 
Of considerable interest is the situation when eo{q) = 0. In this case one obtains 

H [tp - tpis)] =H Atp = il + q)-\ (3.70) 

which generalizes the result of Wainwright [2] for orthogonal spatially homogeneous models exhibiting an 
additional self-similar symmetry. In that case one can choose tp{s) = for ah fluid fines since one has spatial 
homogeneity, thus reducing the above expression to H tp = {1 + q) 

^ Often the notation r is used for proper time which here is denoted by tp. We have followed the conventions of Wainwright 
and Ellis [2] when it comes to the definition of t. Note that this definition differs by a factor of 3 compared to that used by 
Bruni et al [18]. 
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3.6 Discussion 



Above we said that the 6-equations were decoupled, since 6 was the only variable that was not dimensionless. 
However, this is not true if one introduces constants carrying dimension through the matter source or if one 
introduces a cosmological constant (which also carries dimension). Such constants effectively reintroduces 0. 
For example, for a cosmological constant we have that f^A = 3A/0^. This definition leads to the following 
equations: 

BoQa = 2(g+l)nA, (3.71) 
a^^A = 2r„nA. (3.72) 

For a perfect fluid with p = p{p,) as an equation of state it is only p{fi) = (7 — 1) /i that does not lead to 
a reintroduction of 0. For example, for a polytropic equation of state, p(/i) = K fi^ , one obtains f'(fi) = 
3i-')X02(''-i) fi". 

Note also that it is necessary to choose a lapse function Ai (or J\f) and shift M, (or N ) depending only 
on dimensionless variables if one wants to decouple the 0-cquations. This might not be possible if one wants 
to use, for example, comoving local coordinates for the perfect fluid case. Thus dust, for example, leads to 

oc 0. However, in special cases Hke spatially homogeneous models and "silent" models of the Universe 
(to be discussed below) it turns out to be possible to obtain a decoupling of the evolution equation for 0, 
leading to a reduced dynamical problem. Even in the cases where one does not obtain a reduction it might 
still be useful to make a 0-normalization, since this makes it possible to probe the regime where — > ±00. 

The quantity q appearing in the do0-equation can be obtained from Eq. (3.11) (which is the dimensionless 
version of the Raychaudhuri equation). The quantity in the 5Q,0-equation is more problematic. It can be 
obtained from Eq. (3.14) if the matrix ( 2/3 8"^ — E"^ + e"^''' ) is invertible. If one considers the problem 
with two commuting spacelike Killing vector fields, this matrix reduces to a scalar quantity (two components 
of Ta are identically zero if one adapts the frame to the Killing vectors). This is the case considered by 
Hewitt and Wainwright [44]. They claim that this quantity typically is nonzero. Even if the above matrix 
would be degenerate at some point one might still obtain from some of the constraints in the Jacobi 
identities ( these constraints are identically satisfied in the case of two spaceHke KilHng vector fields [44] ) . 
Thus it should probably be possible, generically, to obtain the components of (i.e., those that are nonzero 
even in the presence of Killing vector spacetime symmetries; such symmetries lead to a linear dependence 
between the components Tq, if the symmetry has no timelike component) . 

4 Special Examples 

Perhaps the main application of the above outlined formalism is that it allows one to obtain many interesting 
special cases in a useful unified and compact form just by deleting terms. A number of such cases will be 
discussed in this section. 

4.1 Silent models of the Universe including magnetic Maxwell fields 

In analogy with what is often done in the Newman-Penrose formahsm [15, 48] we can, for example, kine- 
matically speciaHze the preferred timelike congruence u. Here we want to consider nonrotating congruences 
without acceleration, = a-'" = ii". This means that we can impose a synchronous local coordinate system 
with a timelike coordinate affinely parametrizing the congruence. 

Then one can ask the question of what restrictions one needs to impose on the curvature variables in 
order to obtain a decoupled system of evolution equations. This is again in close analogy with the Newman- 
Penrose approach where the restriction to an algebraically special Weyl curvature tensor leads to a decoupling 
of a set of equations which can be solved (see e.g. Refs. [49] and [50] for further references). It turns out 
that in the present case one can obtain a decoupling of a set of evolution equations if one assumes that the 
"magnetic part" of the Weyl curvature tensor vanishes, Hap = 0, and that the matter source defining the 
preferred timelike direction be a nonrotating perfect fluid with dust equation of state, = p = q" = TTap. 
These conditions lead to the generically spatially inhomogeneous "silent" cosmological models of Matarrese, 
Bruni and collaborators [17] - [19]. 
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Here, difTerent from the original models, we consider the possibility of also including a sourcefree magnetic 
Maxwell field, noninteracting with the dust fluid, which is ahgned along one of the eigendirections of the 
"electric part" of the Weyl curvature tensor. Thus we will assume that = Ea = H2 = H3, which, by Eq. 
(2.62), implies that the Poynting vector will be zero. A Maxwell field restricted in this way (together with the 
above mentioned kinematical specialization of the timelike congruence) will have covariantly vanishing spatial 
divergence and spatial rotation [7] . The total stress-energy-momentum tensor of such a matter configuration 
will have the following nonzero contributions: 

H = l^dust + ^ (Hi)'^ , (4.1) 
P = (4.2) 
TTii = -UHif = -27r22 = - 2 7r33 = -i7r+. (4.3) 



Barnes and Rowlingson [51] showed that for irrotational perfect fluid matter with vanishing "magnetic part" 
of the Weyl curvature tensor, a canonical orthonormal frame comoving with u can be chosen such that the 

shear tensor of the fluid congruence, aap, and the "electric part" of the Weyl curvature tensor, Eap, are 
simultaneously diagonahzed. This property foflows from Eq. (2.47), which reduces to the algebraic constraint 



= e''^''apsiE\ + ^n\) , (4.4) 

where we have already taken account of the proposed generalization. The diagonal components of Eqs. (2.41) 
and (2.45) yield that = nn = 7122 = ^133, which in turn implies n% = 0. The off-diagonal components of 
Eqs. (2.27) and (2.44) yield that fi" = 0. Hence, it follows from Eq. (2.11) and = (since a^p and E^p 
can be chosen to be diagonal) that gq, ei, 62 and 63 are HSO (a property which enables one to diagonalize 
the metric tensor), and that the shear eigenframe is Fermi-propagated along u. Thus we have generahzed 
the results of Barnes and Rowlingson [51] to the proposed "silent" configuration with nonzero magnetic 
Maxwell field. We also remark that the momentum conservation equation (2.49) is identically satisfied by 
virtue of the Maxwell field equations. Note also that, since = 1*" = u>", the earlier discussion related to 
inequalities for the values of the deceleration parameter q and the dimensionless product H Atp applies. 

As a result of these specializations, and applying an irreducible decompositions of aap, Eap, iTap and 
*Sap according to Eq. (2.86), the relations presented in Section 2 reduce to the following set: 

The commutators: 

[eo,ei] = -|(0-2a+)ei (4.5) 

[eo,e2] = -^(@ + a++V3a-)e2 (4.6) 

[eo,e3] = -i(0 + a+ -^/3a_)e3 (4.7) 

[ei,e2] = - (02 - risi) ei -I- (ai -I- n23) 62 (4.8) 

[62,63] = - (03 - ni2) 62 -I- (a2 -I- n3i) 63 (4.9) 

[63,61] = - (ai - 7123) 63 -I- (03 -I- ni2) 61 . (4.10) 

Since we have an irrotational geodesic congruence and we demonstrated how a Fermi-transported HSO 
spatial frame {6a} can be chosen, we can now write the orthonormal frame vectors as 

eo = d,^ , 62 = e-^"-^^-^^^- dy , 

where we have used a "Misner parametrization" of the diagonal 3- metric [52]. Of particular interest are the 
commutators involving 60 which give the following relations: 

6o(/3°) = i0, 6o(/3±) = ia±. (4.12) 
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Decoupled subsystem of ordinary difTerential (evolution) equations: 



eo(e) = 


-l@^-2a^-^fidust-k{Hiy'+A 


(4.13) 


eo(cr+) = 


-i(2e -a+)a+-lia-f -E+ + ^iH,f 


(4.14) 


eo((T-) = 


-|(e +a+)a--E_ 


(4.15) 


eo(S+) = 


- i Hdust cr+ - (e + a+) + (T_ B_ + i (e + (T+) 


(4.16) 


eo(£-) = 


- ^ CT- - (e - (T+) + £+ - i 


(4.17) 






(4.18) 


eo(ffi) = 


-|(e + (7+)jfi . 


(4.19) 



Remaining system of evolution equations: 

eo(ai) = -iei(0 + a+)-i(0-2a+)ai (4.20) 

eo(a2) = -ie2(e-ic^+-^c^_)-i(e + c^+ + ^/3(T_)a2 (4.21) 

eo(a3) = -ie3(e-ia+ + ^a_)-i(e + a+-^/3(7_)a3 (4.22) 

eo(n23) = -|(0-2(7+)n23-;)|ei((7_) (4.23) 

eo(n3i) = -i(e + (T+ + ^/3^7_)n3l-ie2((T+-^(T_) (4.24) 

eo(ni2) = -i(e + (T+-^/3^7_)nl2 + ie3((T+ + ^(T_) , (4.25) 

where 



i[(a+)2 + (a_)2] . (4.26) 
Tracefree part and trace of 3-Ricci curvature of spacelike 3-surfaces orthogonal to u (Gaufi equation) : 



*S+ = -i[2ei(ai)-e2(a2)-e3(a3)-4(n23)' + 2(n3i)2+2(ni2)2 
-3(e2 - 2 02) (nai) + 3 (e3 - 203) (7112) ] 

= £;+-i(0 + a+)a+ + i(a_)2 + i(ifi)' (4.27) 
= ^[e2(a2)-e3(a3)-2(n3i)2+2(ni2)2 + 2(ei-2ai)(n23) 

- (e2 - 2 02) (n3i) - (e3 - 2 03) (7112) ] 

= B_ - i (e - 2(T+)(T_ (4.28) 
*R = 2[(2ei-3ai)(ai) + (2e2-3a2)(a2) + (2e3-3a3)(a3) 

- (n23)' - (n3i)' - (ni2)2 ] 

= -lQ^ + 2^ldust + 2<J^ + {H^f+2A. (4.29) 

The constraint equations: 

= ei(0) + (ei -3ai) (ct+) - ^/3n23^T_ (4.30) 

= e2(e)-i(e2-3a2)((T++^/3(7_)-fn3l((7+-^(7_) (4.31) 

= e3(e)-i(e3-3a3)((T+-^/3(T_) + fnl2((T+ + ^(T_) (4.32) 

= (ei-2ai)(n3i) + (e2-2a2)(n23) + ei(a2)-e2(ai) (4.33) 

= (62 - 202) (ni2) + (ea - 203) (risi) + 62(03) - 63(02) (4.34) 

= (63 - 203) (n23) + (ei - 2 oi) (ni2) + 63(01) - 61(03) (4.35) 

= (ei -ai) (ct_) - \/3 n23 0-+ (4.36) 

= (e2-02)(CT+-^CT-) + n3i(CT+ + \/3(T_) (4.37) 

= (e3-03)(f7+ + ^(T_)-nl2((T+-^/3rT_) (4.38) 
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= (ei-3ai)(£;+) + iei(/id„,f)-v^n23£;_ + iai(ffi)2 (4.39) 
= (e2-3a2){E+ + V3E_)-e2{Mdust)+Bn3i{E+-^E_)-l{a-2-n3i}{Hif (4.40) 
= {es-3a3){E+-V3E_)-e3{Mdust)-3ni2{E+ + ^E_)-^{a3+ni2){Hi)^ (4.41) 

(ei - ai) (B_) - V3n23 E+ + ^ 7133 {H^f (4.42) 
= {b2 - a2) {E+ - ^ E_)+ 1131 {E++V2,E_)-\{a2- 1131) {Hif (4.43) 

= {^3-a3){E+ + ^E_)-ni2{E+-V2,E_)-\{a3+ni2){mf (4.44) 

= (ei -2ai)(i/i) (4.45) 

= (e2-a2+n3i)(/?i) (4.46) 

= (e3-a3-ni2)(Hi) . (4.47) 

In order to derive Eqs. (4.39) - (4.44) we made use of the Maxwell equations (4.45) - (4.47). Lesame et al 
[53] recently reported that no new constraints arise from propagating the conditions (4.30) - (4.44) along 
the fluid flow lines in the zero magnetic Maxwell fleld case {Hi = 0). Hence, any Hi = initial data, which 
satisfy the constraints (4.30) - (4.44), evolve according to the evolution equations (4.13) - (4.18) and provide 
a consistent solution to the Einstein field equations with dust matter and = ui" = ii" and Hap = 0. 
Without proof we here assume that this consistency is not being violated by the introduction of the nonzero 
magnetic Maxwell fleld, i?i ^ 0, and leave a rigorous consistency check for a future project. 

It has been conjectured by Berger et al [54] that for irrotational matter flows there might be a link be- 
tween the presence of gravitational radiation in a spacetime and the nonvanishing of the conformal curvature 
of the spacelike 3-surfaces orthogonal to u. Their hypothesis arose from an investigation of the Petrov type D 
spatially inhomogeneous dust spacetimes of Szekeres [55], in which they demonstrated that in those models 
the spacelike 3-surfaces have vanishing conformal curvature. Here we want to reverse the argument. As a 
nonzero "magnetic part" of the Weyl curvature tensor as measured with respect to fundamental observers 
comoving with u is assumed to be a necessary condition for the presence of gravitational radiation in a 
particular spacetime geometry, we are interested in whether the 3-Cotton-York tensor is zero for (general- 
ized) "silent" models and thus further support is given to the Berger et al conjecture. Since the canonical 
orthonormal frame for "silent" models employed here diagonalizes the shear tensor as well as tTq,^ , and only 
off-diagonal terms of the commutation functions Hap are nonzero, it follows directly from Eq. (2.53) that the 
on-diagonal terms of *Ca,8 vanish. However, by use of the (OQ;)-equations (4.30) - (4.32), the i?-constraint 
equations (4.36) - (4.38), and the Maxwell equations (4.45) - (4.47), we obtain from Eq. (2.53) for the 
off-diagonal terms of *Cap 

*Ci = ^,<r- (ei - ai) + n23 [ {a+f - | {cr.f - {Hif ] (4.48) 

*C2 = -^(a+-^a_)(e2-a2)(a_) + v^n3i(a+ + ^a_)a_ + i(a2-n3i)(Ifi)2 (4.49) 

*C3 = (<T+ + ^ (T_) (63 - as) ((T-) - V3ni2 ((7+ - ^ a_) - i (a3 + ni2) {Hif . (4.50) 

We do not see any a priori reason why the righthand sides of these expressions should be zero. Hence, 
against our expectations, for a general representative within the (generalized) class of "silent" models the 
zero "magnetic part" of the Weyl tensor condition does not imply the vanishing of the 3-Cotton-York tensor. 
Nevertheless, the result of Berger et al as regards the dust spacetimes of Szekeres is contained in Eqs. (4.48) 
- (4.50): for = (t_ = n23 ( see Ref. [56] ) and iJi = it follows that = *Ci = *C2 = *C3. 

Dimensionless formulation: 

In terms of the expansion-normalized dimensionless variables introduced in Section 3 the sets of equations 
(4.1) - (4.3), (4.5) - (4.10), (4.14) - (4.19), (4.20) - (4.25), (4.27) - (4.29) and (4.30) - (4.47) become respec- 
tively: 

n = ^dust + ^i'Hi)'' (4.51) 
P = liUif (4.52) 
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(4.53) 
(4.54) 

The commutators: 

(4.55) 
(4.56) 
(4.57) 
(4.58) 
(4.59) 
(4.60) 

On an individual fluid flow line one can choose A'' = ±39, N" = 0. With this choice Eq. (4.12) (with 
which one can replace the above 9o-commutators, if one is so incHned) takes the form 





n+ = 


{Uif 




= 


Q« = n . 


[ao,ai] 


= - n ^0 + {q 


+ 2s+)ai 


[do,d2] 


= -r2do + (q 


-S+ -A/3S_)d2 


[80,83] 


= - rg So + (g 


- s+ + i/3s_)a3 


[81,82] 


= - (ra + 3 ^2 


- 3 A^3l) + (ri + 3 + 3 iV23) ^2 


[82,83] 


= -ir3+SA3 


-37Vi2)a2 + (r2 +3^2 + 3iV3i)a3 


[83,81] 


= -in+SAi 


-3 7V23)a3 + (r3+3A3 + 3iVi2)ai 



ao/3° = 1 , ao/3+ = E± . 

Decoupled subsystem of ordinary differential (evolution) equations: 



80^+ = (g - 1 + S+) E:+ - (S_)2 

aoS_ = (g-l-2E+)S_ -f_ 

8o£+ = (2g-l-3E+)f+-|n 

8o£- = (2g- 1 + 3E+)f_ - In 

8o^dust = (2 (? — 1) fldust 

8oni = {q-i-2^+)ni 

aoOA = 2{q + l) Qa . 



2 ..dust E+ + 3 E_ f_ + I (1 + s+) {niY 



dust 



S_ +3E_f+ - f E_ {-Hif 



(4.61) 



(4.62) 
(4.63) 
(4.64) 
(4.65) 
(4.66) 
(4.67) 
(4.68) 



As before (cf. Eq. (3.7)) we have an evolution equation for 6, which in itself is decoupled from the set 
(4.62) - (4.68): 

8oQ = -il + q)Q , 



and the parameter q is given by 

q = 2[ (E+)2 + (E_)2 ] + i + i (-Hi) 

Remaining system of evolution equations: 

80A1 = qAi + lri-l{8i-ri 
qA2 + ira + i {82 
qA3 + \r3 + l {83 



80A2 
80A3 
80N23 

80N31 

80N12 



■ r2 
r3 



6Ai)E+ 

6.42) (E+ + a/3 E_) 
6A3)(E+ -a/3E_) 
(g + 2E+) 7V23--l=(ai-ri)S_ 

(g - E+ - a/3S_) 7V31 -^82- r2) (S+ - ^ E_) 
(g - E+ + a/3S_) 7Vi2 + I {83 - rs) (S+ + ^ E_) 

Tracefree part and trace of 3-Ricci curvature of spacclike 3-surfaces orthogonal to u (Gaufi equation) : 

<5+ = -H2(ai 



(4.69) 
(4.70) 

(4.71) 
(4.72) 

(4.73) 
(4.74) 

(4.75) 

(4.76) 



n) Ai - {82 - T2) A2 - (83 - r3)A3 - 12 (A^23)' + 6 (7V31)' + 6 (7Vi2)^ 
-3(a2 -r2 -6A2)7V3i +3(a3 -rs -6A3)A^i2 ] 



^+-(l + S+)S+ + (S_)2 + iCHi)2 

^[{82- 7-2) A2 - {83 - r3) A3 -6 (7V3i)2 + 6 (7Vi2)2 + 2 (Si - n - 6 A^) N23 



(4.77) 
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-(a2 

f_-(l-2S+)S_ 
-(251 -2ri - 9Ai)Ai 



r2 -6A2) N31 - {83 - rs - 6 A3) N12 ] 

(253 -2r3- 9^3)^3 



(2^2 -2r2 

^2 



9^2)^2 
2 



— 1 ^dusf 

The constraint equations: 




















+ 3(A^23)' + 3(A^3i)' + 3(A^i2 
[(S+f + (S_f ]-i(?^i)'-i^A . 



2 ^ '31 (S + 



V3 ^- 



) 



(ai + ri)e 
(a2 + r2) e 

(^3 + rs) 

n - (Si - n - 9.4i) S+ + 3\/3 /V23 S 
r2 + I (^2 - 7-2 - 9 .42) (E+ + V3E_) 

rs + 5 (^3 - r3 - 9.43) (E+ - \/3E_) - f 7Vi2 (S+ + ^ S_) 
di - n - 6 Ai) iV3i + (52 - r2 - 6 A2) A^23 + {di - ri) ^2 - (^2 
a2 - r2 - 6 A2) iVi2 + (a3 - r3 - 6 A3) N31 + (82 - r2) A3 - (83 
83 - rg - 6 A3) N23 + (ai - ri - 6 Ai) A^i2 + (33 - rg) A, - (8^ 

81 - n - 3 Ai) E_ - 3^/3iV23 S+ 
a2 - r2 - 3 A2) (E+ - ^ S_) + 3iV3i (S+ + ^/3E_) 
a3 - r3 - 3 A3) (E+ + ^ S_) - 3iVi2 (E+ - ^/3E_) 
Si -2ri -9Ai)f+ + i(ai -2rl)^-3^/3iV23 + f Ai {Tiif 

82 - 2r2 - 9 A2) (£+ + ^/3f_) - {82 - 2r2) + 9iV3i {£+ " 7f ^-) 

f (A2-iV3l)(?^l)' 

a3 - 2r3 - 9 A3) {£+ - ^/3f_) - (33 - 2r3)n - 9A^i2 {£+ + 75^-) 



7-2) Ai 
r3) A2 
n) A3 



f (A3 + iVl2)(?^l)' 

51 -2ri -3Ai)f_ -3^/37V23f+ + ^^^23 (^1)' 
a2 - 2 r2 - 3 A2) {£+ -^£-)+3 N31 {£+ + ^/3 f _) 

a3 - 2 r3 - 3 A3) {£+ + ^£-)-3 N12 {£+ - ^/3 f-) 
di -n -6Ai)-Hi 

52 - r2 -3A2+3iV3i)-Hi 

53 -r3 -3A3 -3Ari2)?^i 

51 -2ri)nA 

52 -2r2)nA 
a3 -2r3)nA . 



1(^2 



N3,){'H,f 



f (A3+iVl2)(?^l)' 



(4.78) 

(4.79) 

(4.80) 
(4.81) 
(4.82) 
(4.83) 
(4.84) 

(4.85) 
(4.86) 
(4.87) 
(4.88) 
(4.89) 
(4.90) 

(4.91) 
(4.92) 

(4.93) 

(4.94) 
(4.95) 
(4.96) 
(4.97) 
(4.98) 
(4.99) 
(4.100) 
(4.101) 
(4.102) 
(4.103) 



The evolution equation (4.68) arises from Eq. (3.71), while the extra constraints (4.80) - (4.82) and (4.101) 
- (4.103) derive respectively from Eqs. (3.8) and (3.72). 

There are a number of ways one can generalize the "silent" models. For example, one can relax the 
conditions on the kinematic quantities and allow rotation. Another possibility is to allow parallel electric 
and magnetic Maxwell fields that are not aligned with the shear or "electric part" eigendirections. The key 
to generalizations seems to be to respect the "silent" conditions: = Hap = q", and to use a geodesic 
congruence, = 0. Perhaps it is not so surprising that, in such cases, one obtains ordinary differential 
equations, since the "silent" conditions are assumed to correspond to the fact that no information is ex- 
changed between different fluid flow lines, motivating the notion "silent" . The above mentioned rotational 
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and electromagnetic generalizations will be discussed further elsewhere [57]. 



4.2 Locally rotationally symmetric models 

LRS symmetry of a fluid spacetime configuration is given when there exists a (normalized) spacelike congru- 
ence e orthogonal to u, covariantly defined by for example a nonzero vorticity vector field, an eigendirection 
of a nonzero degenerate rate of shear tensor field, or a nonvanishing acceleration vector field. As all phys- 
ical measurements have to be invariant under spatial rotations about e and hence identical in all spatial 
directions orthogonal to e, this congruence must be geodesic and shearfree in the local rest 3-spaces or- 
thogonal to u [9, 20]. Given a choice of orthonormal frame with eo = u, ei = e, all covariantly defined 
spacelike vector fields orthogonal to u are consequently colinear with ei , and all covariantly defined spacelike 
symmetric tracefree tensor fields orthogonal to u have coinciding eigenframes with two equal eigenvalues 
(T_ = Ti = T2 = T3 = 0) [9, 20, 58]. Furthermore, one can make a specialization of the orthonormal frame 
such that only the following commutation functions are nonzero [9, 20]: 

0, ill ■= u, cTii = -2(T22 = -2(T33 := -| (^+, 1^1 = ^1 ■= ai := a, 02 = nsi, nu . (4.104) 

The spatial frame derivatives of those variables /, which are proportional to covariantly defined scalar quan- 
tities, have to satisfy the condition = e2(/) = 63(7). 



The commutators: 



eo,ei] = Meo-i(0-2rT+)ei (4.10. 



[eo,e2] = -i(e + rT+)e2 (4.106) 
[60,63] = - I (0 + ^7+) 83 (4.107) 

[61,62] = a 62 (4.108) 
[62,63] = - 2a; 60 + nil 61 + 2 7131 63 (4.109) 
[63,61] = -063 , (4.110) 

where 63(7131) = 0, and rin is proportional to the magnitude of the spatial rotation of the covariantly defined 
spacelike vector field 6. From the LRS spacetime symmetry and a particular use of the freedom of fixing 
the orthonormal frame it follows that 62 and 63 are HSO [9, 20]. In the zero- vorticity LRS case, w = 0, it 
follows that the spatial frame is actually Fermi-propagated along u, and if both = lo = nu, all four frame 
fields are HSO (cf. Eq.(2.11)). 



Table 1: Notational differences. Note that K {"^K) is actually a covariantly defined variable only if a; = = 
nil. 



Stewart and Ellis (1968) [20] van Elst and Ellis (1995) [58] van Elst and Uggla 



a [(0/c) + 2i/3((7/c)]/3 (0-2 (7+)/ 3 

/3 [(0/c)-^/3(a/c)]/3 (0 + ^+)/3 

a —a/2 a 

k —k —nil 

s — 2n3i 

B^K K 2(62-2n3i)(n3i), 
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The evolution equations: 



eo(e) = 


- i 02 + (ei + u - 2 a) (u) - |(ct+)2 + 2 _ 1 + 3p) + A 


(4.111) 




- i (2 - cr+) cr+ - (ei + u + a) (u) + - (E+ - ^ 7r+) 


(4.112) 


69(0;) = 


— 1 (0 + (T+) w + ^ rill ii 


(4.113) 


en fa) 


— ^ (0 + T+) (a + u) — ^ Till ^ ~ \ 


(4.114) 


U \ ' "J- ± / 


- i (0 + 4(T+) nn 

Q ^ 1 ■ '-' -|- J ' " ± _L 


(4.115) 


en(n'l^) = 


— (0 + CT-i-l 


(4.116) 


eo(£+ + i7r+) = 


- i (yii + p) (T+ - (S+ + i 7r+) + i (ei + 2 M + a) (gi) 






-cr+ (£;+ -\t^+) - |niii?+ 


(4.117) 


eo(J?+) = 


- (0 + cr+) lf+ + 1 w gi + 1 mi (£+ - i 7r+) 


(4.118) 


eo(/i) = 


- (/i+p)0-(ei+2w-2a) (gi)-|(T+7r+ 


f4 1 1Ql 


eo(gi) = 


- 1 gi - ei (p) - (/i + p) ii + |(ei + ii - 3 a) (7r+) + | (t+ gi . 


(4.120) 


The constraint equations: 






= (ei-M-2a)(w) 


(4.121) 


= gi-|ei(0)- 


- 1 (ei - 3 a) ((T+) + nil w 


(4.122) 


= (ei - a) (a) + 


1 92 - i (0 + 2 a+)a+ + i (nii)2 + i + i 7r+) - i /i - i A 


(4.123) 


= (ei - 2 a) (rill) - 1 (0 - 2CT+)a; 


(4.124) 


= (ei - a) (risi) 




(4.125) 


= iJ+ + 3 (u + a 


) w + 1 nil cr+ 


(4.126) 


= (ei - 3 a) 


+ i 7r+ ) + i ei (/x) - i (0 + (T+) gi + 3 w J?+ 


(A 1 271 


= (ei-3a)(i?+) + | (/i + p)tj-3tj(£+-i7r+) + fniigi 


(4.128) 


Tracefree part and trace of *Sa(i 


(3-Ricci curvature of spacelike 3-surfaces orthogonal to u when to 


= 0): 


*S+ = 


- ei(a) - [niif + 2 (es - 2 ngi) (rigi) 






+ i 7r+) - i (0 + (7+) cr+ - 3 ^2 


(4.129) 


*R = 


2 (2 ei - 3 a) (a) - i (nii)^ + 4 (ea - 2 risi) (rigi) 






-|02 + 2/i+|(CT+)2 - 6a;2 + 2A . 


(4.130) 



The variables 0, u, a, rin, and correspond directly (up to a constant factor) to covariantly 
defined scalar quantities (cf. Ref. [58]), and therefore, due to the imposed LRS spacetime symmetry, have 
vanishing frame derivatives in the e2- and es-directions. Thus, by appHcation of the commutator relation 
(4.109) to this set of variables, and subsequent substitution from the evolution and constraint equations 
(4.111) - (4.120) and (4.121) - (4.128), one can successively derive consistency conditions which need to be 
satisfied by any solution to the LRS dynamical equations [20]. In this way one obtains from u the condition 

= |(0 + (T+)cj + anii , (4.131) 

which is reproduced when (4.109) acts upon rm; from (0 + a+), using (4.131), the condition 

= cj[(£;+-i7r+) + i(0-cr+)0-|((7+)2-3a;2+3ait+|(nii)2 

+ i(/i + 3p)-A] + |niigi ; (4.132) 

and from a, using (4.131), the condition 

= nii[(£;+ + i7r+) + i(0-CT+)0-|(cr+)2-3a;2 + 3aM+|(nii)2 

- At- A ] -3wgi . (4.133) 
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Combining (4.132) and (4.133) then yields the condition 

= [2w2^i(nn)']gi + [(/i+p)-|7r+]wnii . (4.134) 

This shows that if hypersurface orthogonality properties of the two existing congruences is chosen as a 
classification criterion for consistent LRS spacetime geometries, either cj = or rin = impose strong 
restrictions on the matter field source, since in this case the energy current density is demanded to vanish, 
qi = 0. This is, for example, possible for magnetic Maxwell fields as discussed in the previous subsection. 
No such restriction arises for = oj = nn. Finally, applying (4.109) to E+, we get 

= w[ (/x + p)o-+ - |(e + 2(T+)7r+ - (ei +2M + a)(gi) ] 

+ nn[eiip)-ie + a+)qi] , (4.135) 

while appHed to it reproduces Eq. (4.134). 

Under the assumption that the matter source be a comoving perfect fluid {0 = qi = Tr+) with {fi + p) > 0, 
the consistency condition (4.134) reduces to the simple relation 

= cjnii , (4.136) 

which allows a neat geometric classification of the related LRS spacetime geometries into three distinct 
classes [9, 20]: 

• LRS class 1: CO ^ => = nu = @ = (7-^., eo(/) = . 
eo = u not HSO. 

• LRS class II: = a; = nu . 

All frame vectors HSO. 

• LRS class III: nu ^ ==> = ui = u = a, ei(/) = . 
ei not ESO. 

Solutions within LRS class I are invariant under the transformations of a multiply-transitive G4 with 3-D 
timelike orbits, those within the spatially homogeneous LRS class III under transformations of a multiply- 
transitive G4 with 3-D spacelike orbits. The underlying spacetime symmetry group of LRS class II is 
a multiply-transitive G3 with 2-D spaceHke orbits of either constant positive, zero, or negative Gaufiian 
curvature 

:=2(e2-2n3i)(n3i) . (4.137) 

Contained within the three LRS perfect fiuid classes are a number of prominent cosmological models such as 
the Godel model ( class I: ei(/) = 0, = = (T+ = m = a = nu, = p = A < ), spatially homogeneous 
models (classes II and III: ei(/) =0, = u = ii = a) with Kantowski-Sachs (class II: rin = 0, = 0, 

> ) and Friedmann-Lemaitre-Robertson- Walker ( classes II and III: (t+ = 0, = = ) models as 
special cases, and spherically symmetric Lemaitre-Tolman-Bondi dust models (class II: = a; = a = nu, 
= p = H+, '^K > 0) (sec [9, 20, 58] and references therein). 

Note that Eq. (4.126) implies that LRS class II models have zero "magnetic part" of the Weyl curvature 
tensor. Thus the corresponding dust models are examples of special "silent" models contained in the pre- 
vious subsection. Apart from the Kantowski-Sachs models, further spatially homogeneous specializations 
contained in LRS class II are of LRS Bianchi Type-Ill {^K < 0) and LRS Bianchi Type-I (^K = 0). The 
conformal properties of spacelike 3-surfaces orthogonal to u in LRS models of classes II ( *Cap = ) and III 
( *Cap 7^ ) were investigated by Wainwriglit [59] . Various exact solutions with LRS spacetime symmetries 
have been discussed in Refs. [9], [20] and [48]. 
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Dimensionless formulation: 



In terms of the expansion-normalized dimensionless variables introduced in Section 3 the sest of equations 
(4.105) - (4.110), (4.112) - (4.120), (4.121) - (4.128), (4.129) - (4.130) and (4.131) - (4.135) become respec- 
tively: 



The commutators: 



[do,di] = -{r-3U)do + {q + 2i:+)di 

[80,82] = (g-s+)a2 

[80,83] = (g-E + )S3 

[81,82] = {r + 3A)82 

[82,83] = -6W8o + 3Nii8i+6N3i83 

[83,81] = -{r + 3A)83. 



The evolution equations: 

9oE+ = 

80W -- 

80A -- 

8oNu -- 

80N31 -- 

8o{£+ + ^n+) -- 

80H+ -- 

8ofl -- 

80Q1 -- 

8onA -- 



(g - 1 + E+) S+ - (Si - r + 3 [/ + 3 A) ?7 + 3 VF^ 

(g-l-2S+)VF + f A^nJ7 
qA-U-{U + A)i:+-^NnW-^Qi 
(g-4S+)A^ii 
{q - S+) 7V31 

(2 g - 1 - 3 S+) ^+ + (g + i + i S+) n+ - I (n + P) S+ 

+ ^{8i-2r + 61J + 3A)Qi-^Nnn+ 

(2 g - 1 - 3 E+) 7^+ + I W Qi + f TVn (f+ - i n+) 

(2 g - 1) n - 3 P - (di - 2 r + 6 J7 - 6 A) Qi - 2 S+ n+ 



(f+-in+) 



2(g-i)(3i 
+ 2E+Q1 
2 (g + 1) riA 



[81 - 2r) P - 3 (n + P)U + U8i - 2r + 3U - 9 A) 



As before (cf. Eq. (3.7)) we have an evolution equation for 0, which in itself is decoupled from 
(4.144) - (4.153): 

aoe = -(i + g)0 , 

and 

q=- (di-r + 3U-6A)U + 2 (S+)2 - QW"^ + ^ (fl + 3P) - ■ 
The constraint equations: 



{81 + r) 6 

{81 - r - 3 f/ - 6 .4) W 

Qi + 'Ir - 'lidi - r - 9 A)T.+ + 3 NiiW 

(Si - r - 3 A) A + i - i (1 + 2 S+) S+ + I (7Vn)2 + i (f+ + i n+) - i - i f^A 

{81 - r - 6 A) Nil - 2 {1 - 2T,+)W 

i8i-r-3A)N3i 

n+ + 9{1J + A)W + I Nn^+ 

{8i-2r-9A){£+ + ^ n+) + i (fli - 2r) f2 - f (1 + S+) Qi +9W'H+ 
{81 - 2r - 9 A)n+ + ^ {n + P)W - 9W {£+ - lU+) + ^ Nil Qi 
(Si -2r)0A . 



(4.138) 
(4.139) 
(4.140) 
(4.141) 
(4.142) 
(4.143) 



(4.144) 
(4.145) 
(4.146) 
(4.147) 
(4.148) 

(4.149) 
(4.150) 
(4.151) 

(4.152) 
(4.153) 

the set 
(4.154) 

(4.155) 

(4.156) 
(4.157) 
(4.158) 
(4.159) 
(4.160) 
(4.161) 
(4.162) 
(4.163) 
(4.164) 
(4.165) 
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The evolution equation (4.153) arises from Eq. (3.71), while the extra constraints (4.156) and (4.165) derive 
respectively from Eqs. (3.8) and (3.72). 



Tracefree part and trace of Sap (3-Ricci curvature of spacelike 3-surfaces orthogonal to u when 1^ = 0) : 



S+ = -{di-r)A-3{Nuf + 2{d2-6N3i)Nsi (4.166) 

= (f+ + in+)-(l + E+) (4 1g7) 

/C = -(2di-2r-9A)A+l(N,if -2(d2-6Nsi)Nsi (4.168) 

= 1-n- +9W'^ -flA . (4.169) 

Consistency conditions: 

= lil + T.+)W + ANii (4.170) 
= W[{£+-^U+) + {l-T.+)-2{T.+f-9W'^+9AU+^{Nnf 

+ i(n + 3P)-a\ ] + ! iViiQi (4.171) 
= iVii[ (f+ + in+) + (l-I]+) -2(S+)2-9W^2^9A{7+|(7Vn)2 

-n-nA]-3WQi (4.172) 

= [2W^ + ^iNnf]Qi + [{n + P)-lll+]NiiW (4.173) 
= W[ 3(n + P)E+ - 2(1 + 2S+)n+ - (di - 2r + 6J7 + 3A)Qi ] 

+ 3iVii[(ai-2r)n-(l + S+)Qi] . (4.174) 



4.2.1 LRS class II perfect fluid models 

Although they have been discussed many times in the past (and one will probably continue to do so in 
the future) , we think that in the context of the overall framework presented here it is instructive to take a 
closer look at models within the generically spatially inhomogeneous LRS class II that contain a perfect fluid 
matter source with an equation of state p = p(^i). Using the covariantly defined Gaufiian 2-curvature of the 
symmetry group orbits, "^K, in favor of E+ as a dynamical variable, the sets of equations (4.105) - (4.110), 



(4.111) - (4.120), and (4.121) - (4.128) specialize to: 

[eo,ei] = ueo - i (6 - 2(T+)ei (4.175) 

[eo,e2] = -i(0 + a+)e2 (4.176) 

[eo,e3] = -|(e + a+)e3 (4.177) 

[ei,e2] = a 62 (4.178) 

[62,63] = 2 713163 (4.179) 

[83,81] = -aeg , (4.180) 

eo(e) = -ie2 + (ei+u-2a)(u)-|(cr+)2-i(;i + 3p) + A (4.181) 

6o((7+) = -\Q'^ -{Q-\a+)a+-{ei+u + a){u) + la^ -\'^K+\n+\K (4.182) 

eo(«) = ei(f 0) + (f -i)0« + |a+« (4.183) 

eo(a) = -\{Q + a+){a + u) (4.184) 

eo{^K) = -|(e + cr+)2X (4.185) 

eo(/i) = -(/i + p)e, (4.186) 

(Note that we have added a ii evolution equation since we are considering a nontilted perfect fluid.) and 

= ei(0) + (ei -3a) (cr+) (4.187) 

= (ei -|a)(a) + i02-i(CT+)2 + i2;f-|/i-|A (4.188) 

= {ei-2a)CK) (4.189) 

= ei{p) + {n+p)u . (4.190) 
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E+ and are related by 

= f 2^ + i (e + a+f - fa^-i/i-iA. (4.191) 

Now taking the slicing point of view we can choose local coordinates t,x,y,z and express the orthonormal 
frame vectors according to Eq. (2.76) as 

eo = (dt - d,) , 62 = dy, 

ei=X-iS, , B^ = {YZ)-^d,, ^^■'''^> 

where iV = A^{t,x), A^^ = Nx{t,x), X = X{t,x), Y = Y{t,x), and Z = Z{y). Inserting these expressions 
into the commutation relations (4.175) - (4.180) we obtain 

u = N-'^eiiN) (4.193) 

eoiX) = l{e-2a+)X + X''N-^eiiN^) (4.194) 

eo(F) = i(e + CT+)y (4.195) 

= (ei+a)(y) (4.196) 

= (e2 + 2n3i)(Z) . (4.197) 

Then, demanding that '^K be constant for t = const = x, we get from Eq. (4.137) the condition 

= Z^yy + CiZ , Ci= const , (4.198) 

from which it follows that "^K = Ci / Y"^ . If Ci > the group orbits are spherically symmetric 2-surfaces, if 
Ci = they are flat 2-planes, and if Ci < the 2-surfaces have hyperbolic geometry. 

How much information about the metric docs the system of equations (4.181) - (4.190) contain? Through 
the commutators we flnd that the metric variable Y occurs implicitly through in the spherically and hy- 
perbolically symmetric cases. However, there is no information about X, which is impHcitly appearing in 
these equations through ei . Thus to obtain a complete system one has to include the commutator equation 
corresponding to Eq. (4.194). It is also worth noting that the time and space frame derivatives of Y are 
encoded through the kinematic quantities (by the combination (0 + (t+) ) and the variable a. This is in 
stark contrast to the metric variable X, where the only information we can obtain is that its eo-derivative is 
determined by the combination (6 — 2 (t_|_), if we have specified a shift vector (of course the shift also occurs 
implicitly in gq in the above equations). The fact that we do not have a variable corresponding to an ei- 
derivative of X is completely analogous to the nonexistence of an eo-derivative of the lapse; the nonexistence 
of an ei -derivative of X is associated with the freedom of making a coordinate reparametrization x = x{x). 

In terms of dimensionless variables the LRS class II perfect fluid models are described by the sets of equations: 



[do,di] = -(r-3U)do + (q + 2i:+)d, 
[80,82] = {q-^+)82 

[80,83] = (g-S+)a3 

[31,^2] = {r + 3A)82 

[82,83] = 6N3183 

[83,81] = -{r + 3A)d3, 

aoS+ = (g-2 + iS+)S+-i-(ai-r + 3?7 + 3A)J7+f A2_|2^^if^^i^^ 

9nU = (g + 3f +2S+)J7 + (ai-r)(f ) 

80A = qA-U - (il + A)'E+ 

ao% = 2 (g - S+) % 

8ofl = {2q-l)n-3P 

8onA = 2 (g + 1) , 



4.199) 
4.200) 
4.201) 
4.202) 
4.203) 
4.204) 

4.205) 
4.206) 

4.207) 
4.208) 
4.209) 
4.210) 
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where 
and 



q = - (di -r + 3U -6A)U + 2 (S+)2 + i (fi + 3P) _ fi^ , 



(4.211) 



= r - (di - r - 9 A) S+ 

= (ai-r-|A)A + i-i(S+f + i% 

= {di-2r-6A)^IC 

= {di-2r)P + 3{n + P)U 

= (Si - 2 r) r^A ■ 



(4.212) 
(4.213) 
(4.214) 
(4.215) 
(4.216) 



The variables £+ and % := 3'^K / 0^ are related by 



|% + i(l + S+)2-|A2-in-inA. 



(4.217) 



Let us now consider the special case where p(fi) = (7 — 1) /U, P{^) = (7 — 1) ^, so that is not "rein- 
troduced". Equation (4.211) defines q while one can use Eq. (4.213) to solve for f!. One can choose Eq. 
(4.212) to solve for r. This is indeed the equation Hewitt and Wainwright used in the case of two commuting 
spacelike Killing vector fields (which overlaps with the present case when = 0) [44]. However, in the 
present case one can also choose (4.214) in order to solve for r. Once one has solved for r the remaining 
equation will constitute a constraint (in contrast to the case considered by Hewitt and Wainwright which 
did not contain such a constraint [44]). If one includes a cosmological constant, Eq. (4.216) implies that the 
possibilities of solving for r increase even further. 

4.3 Spatially homogeneous models 

When dealing with spatially homogeneous models it is convenient to choose eo as the unit normal with 
respect to the spacelike 3-surfaces of homogeneity. Since one usually focuses on spatially homogeneous 
expanding models, in this section we will go directly to the dimensionless formulation. The condition of 
adapting to the spatially homogeneous 3-surfaces leads to the restrictions: = 11°' = = Ta = daF; F 
denoting the dimensionless form of any geometrically defined variable. (Note that, since = U" = W", the 
earlier discussion related to inequalities for the values of the deceleration parameter q and the dimensionless 
product H Atp applies.) Spatially homogeneous models have been reviewed many times before, see, for 
example, Refs. [21], [22], and [2]. Since we have adapted to the spatially homogeneous 3-surfaces it is 
natural to take the slicing point of view. This leads to the relations 



where e, are the vector fields associated with the spatial homogeneity symmetry group. Thus ^''jk are just 
the structure constants describing the various Bianchi groups (the Kantowski-Sachs models with multiply- 
transitive spacetime isometry group have already been treated in the previous LRS section). Spatial homo- 
geneity requires Af = M{t), while the shift vector can be set to zero (however, a full use of the so-called 
automorphism group requires a nonzero spatially dependent shift vector satisfying certain requirements; for 
further discussions see Refs. [60] and [61] ). In this case one can choose A/" = 3, leading to a time variable 
T = ln(£/£o) (recall Eq. (1.8)). Thus the operator do acting on a function F{t) reduces to dF/dr. These 
restrictions lead to: 

The commutators 



do = 



■I 1 



(4.218) 



[do, da] 



3 [ig<5^„-E^„+e^„^i?^] dp 

3 [2A^aPp]+ea,06N''^] d^ 



(4.219) 
(4.220) 
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The field equations 

Decoupled equation 



doQ = -{l + q)Q, (4.221) 
q = 2S2 + i(n + 3P)-nA ■ (4.222) 



Remaining field equations 



SoE"^ = (g-2)E"^+n"^-<S"^ + 6eT"'("iZ^S^^5 (4.223) 
n = 1 - - /C - f^A (4.224) 
= Q" -9Api:"'^ -3e"'^"' Npsi:% , (4.225) 



where 



Sap = SBa0-SapB\+Qe'\^A\^\Np)s (4.226) 

K. = 9AaA" + ^B"a (4.227) 

Bap = 2Na-,N^p-N^^Nap (4.228) 

The Jacobi identities 



doA" = qA" -3Ap Y."^ - 3 e"'^'' A13 i?., (4.229) 

doN"'^ = qN"^ + 6T,'^°'^N'^^"' -6e"'^'-"N%Rs (4.230) 

= ApN"^ (4.231) 

The "electric" and "magnetic parts" of the Weyl tensor 



^aP + ^ ^aP — ^aP - 3 ^^P + | SaP + SaP (4.232) 

■Hap = lN\i:ap-9N\„Yp^^ + 3dapN^s^''' -Se'^'iaA^^i^p^s (4.233) 
The Bianchi identities for the Weyl tensor 



ao(f"^ + in"^) = (2(7-1)^"^ + (g + i)n"^- f (o + p)s"^-|A("Q^) 

+ 9 E% (f '^^'^ - i n^)'^) + f N\ n"'^ - 9 7v("^ -H^)-^ 

+ [ i Q'' - 3 (^-^-^ - i n-"*) + 3 iV^a -H'^'' ] 

- 3 e^^(« [ A, T^'^),- - 2 i?^ (£'^\ + i n'^)^) - 5 (9^ ] (4.234) 

aoi?"^ = (2g-i)?^"^ + 9s(%-H^)T-|Ar'^^(i:"'^-in"^)+ 9 iv^"^ (f - i n^)^) 

- 3 S-^^ [ Ts^s + N^s {£^^ - \ W^) ] 

+ 3 [ A^) [£^\ - \ n^),) + \ S^)^ Q, + 2 i?^) T^''), ] (4.235) 

Note that the constraints that arise in tiie Bianciii identities do not give any new information. They are just 
combinations of the constraints obtained in the field equations and the Jacobi identities. 



34 



The Bianchi identities for the source terms 



doil = (2g-l)n-3P + 6AaQ"-3E„^n"^ (4.236) 
doQ" = 2iq-l)Q" + 9ApIl"^ -3 + 3 e«^^ [RpQ-, + Nps ] ■ (4.237) 

The dimensionless form of the 3-Cotton-York tensor for spatially homogeneous models is 

Cap = - (So - 2 g + 2) Hap + 9 E^'(„ Hp)^ - 27 Y.% + f TV^ ^'p + 6 iV„^ E^ 

- 3 (5„/3 [ - 3 A^^a S*, E% + E^ - 3 A^^a ] 

- 3 [ (3 E^^) E,, + n^),) - i E^)^ Qs - 2 J?,^, ?^^), ] . (4.238) 

The conformal properties of spacelike 3-surfaces in (orthogonal) spatially homogeneous cosmological models 
were investigated by Wainwright [59]. 

Note that the class B parameter h gives rise to a constraint through its definition: 

AaAp = ^hea^,epsrN''^ N"- ^ = \hN , N := {N\f - N^pN"!^ . (4.239) 

Expressed in expansion-normalized variables the Bianchi-Schiicking-Behr classification ( see e.g. Ref. [22] ) 
takes the form given in Tab. 2 (the variables, and therefore also the present form of the classification, breaks 
down when 6 = 0. In that case one has to use other variables, such as for example the original dimensional 
ones) . It is interesting to note that the Type-I models ( which are characterized by = Nap = A" ) are 
examples of "silent" models (see Eq. (4.233)). 



Table 2: Classification of the Bianchi groups into classes A and B and group types I-IX. Bianchi Type-Ill 
belongs to Type-VI^ with h = —1. 



Group Class 




Group Type 




Class A Aa = 


det {N"l^) ^ 


IX: sgn (6) det (iV"^) > 


Vni: sgn (6) det (iV"^) < 




det (iV«^) = 


VIIo: N>0 


VIo: N <0 






II: TV = , 7V"„ ^ 


I: N = 0, N% = 


Class B Aa^O 


det (iV«^) = 


Vllft: N>0 


Yh: N <0 






IV: N = 0,N"aj^O 


Y: N = 0,N'^a = 



Note that the above equations are ordinary differential equations. Furthermore, at first sight there seems 
to be no need for the commutators, since the only differential operator appearing in them (flo) reduces to an 
ordinary derivative (d/dr), which is completely known. This is a consequence of the spatial symmetry. In 
more general circumstances one does not have control of the information hidden in the differential operators. 
Thus in such cases one is forced to consider the commutator equations. However, even in the spatially 
homogeneous case one needs to consider the commutators if one wants to construct the metric. Nevertheless, 
in this case it turns out that one can do so algebraically without solving any differential equations, by choosing 
a shift vector adapted to the automorphism group [61]. The most common frame choice is one restricting 
Nap and A". Usually, with the exception of some special class B models, one chooses Nap to be diagonal. 
Various exact solutions of spatially homogeneous (fiuid) spacetime geometries have been discussed in, for 
example, Refs. [21], [22], [48] and [62]. 
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4.3.1 Perfect fluids 



Here we will consider a perfect fluid with equation of state p{jl) = (7 — 1) fl. Since we want to assume 
that the fluid is tilted with respect to the spatially homogeneous 3-surfaces, this leads to the dimensionless 
expressions given in Eq. (3.37). Then the equations which follow for fl and v" are 



^0^" = r7T-7I-TT::JT [5(37-4)(l-«^) + (2-7)(S^7«^«^) 



don = fiG-i [ 2Gg- (37- 2) - (2 -7)^2 _3-yE„^w"w^ +67A„w" ] (4.240) 

7[ 1-F^1)«2 ] 

+ [(3-27) + (7 -l)v'] (An v^) ] - 3 S"^ + 3 e«^^ [Rp + N^p vs) 
-Sv^A". (4.241) 

It is convenient to use Eq. (4.224) to solve for fJ. However, it is often of advantage not to use the constraint 
(4.225) to solve for v" globally, since the constraint surface is quite compHcated. Instead one can use Eq. 
(4.225) to solve for some variables (not necessarily v") locally. For examples on how to deal with constraints 
of the type given in equation (4.225), see Refs. [63], [64] and [65]. It is usually convenient to make variable 
transformations to other dimensionless variables when dealing with specific problems. Examples of this can 
be found in Ref. [66], which addressed class B models, and Ref. [46], which dealt with LRS Type-IX models. 

4.3.2 Maxwell vacuum flelds 

The dimensionless form of the contributions to the Maxwell stress-energy-momentum tensor were given in 
Eqs. (3.41) - (3.43), while the sourcefree Maxwell equations for spatially homogeneous fields are given by 



dnS" = {q-l)S" +3T,°'p£'^ -3N"p'H'^ -Se"'^^ [Ap-H-y-RpE-y] (4.242) 

ao-H" = {q-l)'H°' + 3T,''p'H'^ +3N"p£'^ + 3e'"^"' [ApE^ + Rp-H^] (4.243) 

= Ac£" (4.244) 

= Ac-H" . (4.245) 



5 Concluding Remarks 

It has been shown how, from the general 1+3 orthonormal frame equations presented in sections 2 and 3, 
one can extract equations suitable for further investigations for a large set of different problems. In a series 
of future papers the authors will, together with M. Bruni, present a dynamical systems analysis for rotating 
"silent" models and "silent" models including electromagnetic fields [57]. 

The dimensionless formulation presented in section 3 is adapted to the existence of homothetic Killing 
vector fields, since the result of acting with the associated symmetry transformation on any dimensionless 
quantity yields zero. It is believed that self-similar solutions act as key building blocks when it comes 
to the understanding of more general models not exhibiting such a symmetry. For example, this is true 
in spatially homogeneous cosmology and for models with two spacelike commuting Killing vector fields, 
where they occur as equilibrium states in the reduced dimensionless phase space (see e.g. Ref. [2]). Thus 
dimensionless formulations should provide a proper setting for probing the role of self-similar solutions in 
General Relativity. 

The 1+3 formulation which has been presented here is motivated by the existence of a preferred timelike 
vector field u. Usually this field is supplied by the stress-energy-momentum tensor. When it comes to vacuum 
spacetime geometries it is perhaps more natural to use formahsms based on 2 + 2 splittings rather than 1+3 
or 3 + 1 splittings. In a 2 + 2 splitting formalism one can choose to adapt to null congruences. This has the 
advantage that the formalism then is naturally adapted to the algebraic properties of the Weyl curvature 
tensor. Moreover, such approaches should also be useful for problems involving gravitational radiation. In 
this context one can mention the Newman-Penrose formalism and the Geroch-Held-Penrose formalism ( see 
e.g. Refs. [67], [14] ). Another 2 + 2 formalism worth mentioning is the covariant 2 + 2 formalism developed 
by d'Inverno and others [68], which now frequently finds application in the literature. Here the spacetime 
manifold is foliated by two families of spacelike 2-surfaces, which are Lie-dragged along two congruences of 
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arbitrarily timelike, null or spacelike character. Its advantages lie in the fact that it can explicitly identify 
the degrees of freedom of a vacuum gravitational field as the entries of a 2-metric, and that initial data is 
unconstrained. 

In general 2 + 2 approaches seem to be preferred when it comes to vacuum spacetimes. However, it should 
be pointed out that there exist nontrivial vacuum and electromagnetic invariant submanifolds in the "silent" 
cases, particularly if one allows the congruence u to have rotation. Moreover, if (t_ and are nonzero 
one has Petrov type I. Since most exact vacuum and electromagnetic solutions are algebraically special, the 
models arising in the 1+3 context, which lead to a decoupled system of ordinary differential equations, may 
provide an alternative to the traditional 2 + 2 approaches when it comes to finding out properties of models 
that are of Petrov type I. 

In the context of "silent" models one may investigate the size of the solution space for a given proposed 
configuration. Some insight in this respect can be gained from a systematic study of the integrability 
conditions underlying the particular setting under scrutiny, that is, various consistency conditions might arise 
from propagating the constraint equations along the preferred timelike congruence u. Another interesting 
possibility of probing the solution space is provided by an application of the equivalence problem approach 
[38, 39]. Further progress could arise from numerical calculations aimed at revealing the content of the 
related constraint equations, which, to our knowledge, to date have not been pursued. Finally, it should 
be pointed out that "silent" models contain a fairly large set of known exact solutions, although no new 
solution has been discovered within this class so far. However, it seems unhkely that there should not 
exist any further solution for these models, apart from the known exact ones. Even if this was so, the 
given geometric formulation of the class of "silent" models (and beyond) presents a remarkably attractive 
framework that nicely unifies otherwise scattered results. 
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A Appendix 

When translating the equations of the 1+3 splitting approach to General Relativity from covariant into 
orthonormal frame form the fohowing relations vahd for the derivative terms occurring are very useful: 



A.l Vector derivatives 



ui" = 0. 



h^^uPVpV =^ eo{V")-e'"^^npV-, . (A.l) 



h'^^V^V =^ (e„-2a„)(y") . (A.2) 
^^,upa (y^y^) ^ _a^^v^)- n"p . (A.3) 
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A. 2 Tensor derivatives 

h''ph''„u''VrAP'' => eo{A"'^)-2e'^^^"n^A'^\ . (A. 5) 

h% hW.AP" =^ (e;3 - 3 ap) {A"^) - e"^'' n^s A'^ . (A.6) 

h^^p h''\ T^P^^^iVrA"^) ux =^ e^^^" (e^ - a^) (Af^^) -3n% A^^^"' + i n\ yl"^ + 5"^ n''s A\ . (A.7) 

References 

[1] J. Ehlers, "Beitrage zur relativistischen Mechanik kontinuierlicher Medien", Akad. Wiss. Lit. Mainz, 
Abhandl. Math.-Nat. Kl 11, (1961). 

Translation: J. Ehlers, "Contributions to the Relativistic Mechanics of Continuous Media" Gen. Rel. 
Grav. 25, 1225 (1993). 

[2] J. Wainwright and G. F. R. Elhs (Eds.), Dynamical Systems in Cosmology (Cambridge University Press, 
Cambridge, 1996). To appear. 

[3] C. W. Misner, K. S. Thorne and J. A. Wlicclcr, Gravitation (Freeman & Co., New York, 1973). 

[4] G. F. R. Ellis, "Cosmological Models from a Covariant Viewpoint", University of Cape Town Preprint 
(1996). Talk given at ICGC95. To appear in Conference Proceedings. Ed. Padmanathan et al. 

[5] S. W. Hawking, "Perturbations of an Expanding Universe", Astrophys. J. 145, 544 (1966). 

[6] G. F. R. Ellis, "Relativistic Cosmology", General Relativity and Cosmology Proceedings of the XLVII 
Enrico Fermi Summer School, Ed. R. K. Sachs (Academic Press, New York, 1971). 

[7] G. F. R. Ellis, "Relativistic Cosmology", Cargese Lectures in Physics Vol. 6, Ed. E. Schatzman (Gordon 
and Breach, New York, 1973). 

[8] F. A. E. Pirani, "On the Physical Significance of the Riemann Tensor", Acta Phys. Polon. 15, 389 
(1956). 

[9] G. F. R. Ellis, "Dynamics of Pressure- Free Matter in General Relativity", J. Math. Phys. 8, 1171 (1967). 

[10] M. A. H. MacCallum, "Cosmological Models from a Geometric Point of View", Cargese Lectures in 
Physics Vol. 6, Ed. E. Schatzman (Gordon and Breach, New York, 1973). 

[11] M. A. H. MacCallum, private notes. 

[12] R. M. Wald, General Relativity (University of Chicago Press, Chicago, 1984). 

[13] F. de Fehce and C. J. S. Clarke, Relativity on Curved Manifolds (Cambridge University Press, Cam- 
bridge, 1990). 

[14] S. B. Edgar, "The Structure of Tetrad Formalisms in General Relativity: The General Case", Gen. Rel. 
Grav. 12, 347 (1980). 

[15] E. Newman and R. Penrose, "An Approach to Gravitational Radiation by a Method of Spin Coeffi- 
cients", J. Math. Phys. 3, 566 (1962). 

[16] R. T. Jantzen, P. Carini and D. Bini, "The Many Faces of Gravitoelectromagnetism" , Ann. Phys. 215, 
1 (1992). 



38 



[17] S. Matarrese, 0. Pantano and D. Saez, "General Relativistic Dynamics of Irrotational Dust: Cosmolog- 
ical Implications", Phys. Rev. Lett. 72, 320 (1994). 

[18] M. Bruni, S. Matarrese and O. Pantano, "Dynamics of Silent Universes", Astrophys. J. 445, 958 (1995). 

[19] M. Bruni, S. Matarrese and 0. Pantano, "A Local View of the Observable Universe", Phys. Rev. Lett. 
74, 1916 (1995). 

[20] J. M. Stewart and G. F. R. Ellis, "Solutions of Einstein's Equations for a Fluid which Exhibits Local 
Rotational Symmetry", J. Math. Phys. 9, 1072 (1968). 

[21] G. F. R. Ellis and M. A. H. MacCallum, "A Class of Homogeneous Cosmological Models", Commun. 
Math. Phys. 12, 108 (1969). 

[22] A. R. King and G. F. R. Ellis, "Tilted Homogeneous Cosmological Models", Commun. Math. Phys. 31, 
209 (1973). 

[23] R. T. Jantzen, P. Carini and D. Bini, Understanding Spacetime Splittings and Their Relationships, in 
preparation. 

[24] P. Szekeres, "The Gravitational Compass", J. Maths. Phys. 6, 1387 (1965). 

[25] J. L. Synge and A. Schild, Tensor Calculus (University of Toronto Press, Toronto, 1949). 
Reprinted: (Dover Publ., New York, 1978). 

[26] W. Kundt and M. Triimper, "Beitragc zur Theorie der Gravitations-Strahlungsfelder" , Akad. Wiss. Lit. 
Mainz, Abhandl. Math.-Nat. Kl. 12, (1962). 

[27] E. Cotton, Ann. Fac. Sci. Toulouse (II) 1, 385 (1899). 

[28] J. W. York, Jr., "Gravitational Degrees of Freedom and the Initial- Value Problem", Phys. Rev. Lett. 
26, 1656 (1971). 

[29] H. Stephani, Allgemeine Relativitdtstheorie (4. Aufl.) (Dt. Verlag d. Wissenschaften, BerHn, 1991). 

Translation: H. Stephani, General Relativity (Cambridge University Press, Cambridge, 2nd Edn. 1990). 

[30] S. Boersma and T. Dray, "Slicing, Threading and Parametric Manifolds", Gen. Rel. Grav. 27, 319 (195). 

[31] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Pergamon Press, Oxford, 4th Edn. 
1975). 

[32] A. H. Taub, "Stability of Fluid Motions and Variational Principles", Proceedings of the 1967 Colloque 
on "Fluids et Champ Gravitationel en Relativite Generale," No. 170, (Centre National de la Recherche 
Scientifiquc, Paris), 57 (1969). 1 (1992). 

[33] J. W. York, Jr., "Kinematics and Dynamics of General Relativity", Sources of Gravitational Radiation, 
Procceedings of the Battelle Seattle Workshop, Ed. L. L. Smarr (Cambridge University Press, Cambridge, 
1979). 

[34] W. Israel, "Nonstationary Irreversible Thermodynamics: A Causal Relativistic Theory", Ann. Phys. 
(N.Y.) 100, 310 (1976). 

[35] W. Israel and J. M. Stewart, "Transient Relativistic Thermodynamics and Kinetic Theory", Ann. Phys. 
(N.Y.) 118, 341 (1979). 

[36] A. Papapetrou, "Quelques Remarques sur le Formalisme de Newman-Penrose" , C. R. Acad. Sci. (Paris) 
A 272, 1537 (1971). 

[37] A. Papapetrou, "Les Relations Identiques entre les Equations du Formalisme de Newman-Penrose" , C. 
R. Acad. Sci. (Paris) A 272, 1613 (1971). 



39 



[38] E. Cartan, Legons sur la Geometrie des Espaces de Riemann (Gauthier-Villars, Paris, 2nd Edn. 1946). 

[39] A. Karlhede, "A Review of the Geometrical Equivalence of Metrics in General Relativity" , Gen. Rel. 
Grav. 12 (1980), 693. 

[40] M. Bradley and A. Karlhede, "On the Curvature Description of Gravitational Fields", Class. Quantum 
Grav. 7 (1990), 449. 

[41] H. Bondi, "The Contraction of Gravitating Spheres", Proc. Roy. Soc. London A 281, 39 (1964). 

[42] B. K. Harrison, K. S. Thorne, M. Wakano and J. A. Wheeler, Gravitation Theory and Gravitational 
Collapse (The University of Chicago Press, Chicago and London, 1965). 

[43] C. B. Collins, "More Qualitative Cosmology", Commun. Math. Phys. 23 (1971), 137. 

[44] C. G. Hewitt and J. Wainwright, "Orthogonally Transitive G2 Cosmologies", Class. Quantum Grav. 7, 
2295 (1990). 

[45] K. Rosquist and R. T. Jantzen, "Unified Regularization of Bianchi Cosmology", Phys. Rep. 166, 89 
(1988). 

[46] C. Uggla and H. von Zur-Miihlen, "Compactified and Reduced Dynamics for Locally Rotationally 
Symmetric Bianchi Type IX Perfect Fluid Models", Class. Quantum Grav. 7, 1365 (1990). 

[47] D. M. Eardley, "Self-Similar Spacetimes: Geometry and Dynamics", Commun. Math. Phys. 37, 287 
(1974). 

[48] D. Kramer, H. Stephani, M. A. H. MacCallum and E. Herlt, Exact Solutions of Einstein's Field Equa- 
tions (VEB Dt. Verlag d. Wissenschaften, Berlin, 1980). 

[49] J. Stewart, Advanced General Relativity (Cambridge University Press, Cambridge, 1990). 

[50] G. Ludwig and S. B. Edgar, "Integration in the GHP Formahsm L A Coordinate Approach with 
Apphcations to Twisting Type N Spaces". To appear in Gen. Rel. Grav. (1996). 

[51] A. Barnes and R. R. Rowlingson, "Irrotational Perfect Fluids with a Purely Electric Weyl Tensor", 
Class. Quantum, Grav. 6, 949 (1989). 

[52] C. W. Misner, "The Isotropy of the Universe", Astrophys. .J. 151, 431 (1968). 

[53] W. M. Lesame, P. K. S. Dunsby and G. F. R. Ellis, "Integrability Conditions for Irrotational Dust with 
a Purely Electric Weyl Tensor: A Tetrad Analysis", Phys. Rev. D 52, 3406 (1995). 

[54] B. K. Berger, D. M. Eardley and D. W. Olson, "Note on the Spacetimes of Szekeres", Phys. Rev. D 16, 
3086 (1977). 

[55] P. Szekeres, "A Class of Inhomogeneous Cosmological Models", Commun. Math. Phys. 41, 56 (1975). 

[56] D. A. Szafron and C. B. Collins, "A New Approach to Inhomogeneous Cosmologies: Intrinsic Symme- 
tries. II. Conformally Flat Slices and an Invariant Classification", J. Math. Phys. 20, 2354 (1979). 

[57] M. Bruni, H. van Elst and C. Uggla, in preparation. 

[58] H. van Elst and G. F. R. Elhs, "The Covariant Approach to LRS Perfect Fluid Spacetime Geometries", 
gr-qc/9510044, QMW/UCT Preprint, (1995). To appear in Class. Quantum Grav.. 

[59] J. Wainwright, "A Classification Scheme for Non-Rotating Inhomogeneous Cosmologies", J. Phys. A: 
Math. Gen. 12, 2015 (1979). 



40 



[60] R. T. Jantzen, "Spatially Homogeneous Dynamics: A Unified Picture", in Proc. Int. Sch. Phys. "E. 
Fermi" Course LXXXVI on "Gamov Cosmology" (R. Ruffini, F. Melchiorri, Eds.) (North Holland, 
Amsterdam, 1987), and in Cosmology of the Early Universe (R. Ruffini, L.Z. Fang, Eds.) (World 
Scientific, Singapore, 1984). 

[61] C. Uggla and R. T. Jantzen, in preparation. 

[62] C. Uggla, R. T. Jantzen and K. Rosquist, "Exact Hypersurface-Homogeneous Solutions in Cosmology 
and Astrophysics", Phys. Rev. D 51, 5522 (1995). 

[63] C. G. Hewitt and J. Wainwright, "Dynamical Systems Approach to Tilted Bianchi Cosmologies: Irro- 
tational Models of Type V", Phys. Rev. D 46, 4242 (1992). 

[64] U. Nilsson and C. Uggla, "Spatially Self-Similar Locally Rotationally Symmetric Perfect Fluid Models", 
gr-qc/9511064. Submitted to Class. Quantum Grav.. 

[65] U. Nilsson and C. Uggla, "Self-Similar Spherically Symmetric Perfect Fluid Models", in preparation. 

[66] C. G. Hewitt and J. Wainwright, "A Dynamical Systems Approach to Bianchi Cosmologies: Orthogonal 
Models of Class B", Class. Quantum Grav. 10, 99 (1993). 

[67] R. Penrose and W. Rindler, Spinors and Space-Time, Vol. I (Cambridge University Press, Cambridge, 
1984). 

[68] R. A. d'Inverno and J. Smallwood, "Covariant 2-1-2 Formulation of the Initial- Value Problem in General 
Relativity", Phys. Rev. D 22, 1233 (1980). 

* * * 



41 



